Weibull in BRMs

Julin Maloof

Weibull in brms

Data set:

Use all dates for the Weibull modeling

dat <- read_csv("../output/height_data_clean.csv")

dat <- dat %>%
  mutate(day=as.numeric(survey_date-min(survey_date)+30),
         mf=factor(mf))

datbcd <- dat %>% filter(str_detect(pop, "BH|CC|DPR"))

Non-linear model specification in brms

Our formula specification is a bit more complex, so we use the function brmsformula to define the non-linear formula.

I am also going to rename a couple of the parameters to be more clear

f1 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmax + Hmin + k + delta ~ 1,
                  nl=TRUE)
  • The first line should look familiar (Weibull formula)
    • I am dividing k by 100 just to get it on a more reasonable scale
  • The second line is new, and it defines how each parameter should be fit with respect to the predictors in our data set.
  • The notation Hmax + Hmin + k + delta ~ 1 means that we want a single parameter value (Intercept) for each of the four parameters
  • nl=TRUE tells brms that this is a non-linear formula.

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)

Set the priors

What priors should we use?

prior1 <- c(set_prior("normal(5,3)", nlpar="Hmin"),
            set_prior("normal(60,15)",  nlpar="Hmax"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta")
)

Why gamma?

  • Want to constrain k and delta to be > 0.
  • Setting lower bounds causes problems
  • Gamma distribution is positive
tibble(x=seq(0,10,.1),
       gamma_2_2=dgamma(x, 2, 2),
       gamma_20_3=dgamma(x, 20, 3)) %>%
  pivot_longer(-x, values_to = "density") %>%
  ggplot(aes(x=x,y=density, color=name)) +
  geom_line()

Fit the model

fit1 <- brm(formula=f1, data=datbcd, prior=prior1, cores = 4)
Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1

Get the summary

summary(fit1)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmax ~ 1
         Hmin ~ 1
         k ~ 1
         delta ~ 1
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmax_Intercept     65.45      1.49    62.76    68.57 1.00     2266     2355
Hmin_Intercept      4.23      0.54     3.16     5.25 1.00     2852     2893
k_Intercept         0.67      0.01     0.64     0.69 1.00     2227     2273
delta_Intercept     6.09      0.48     5.20     7.08 1.00     2195     2370

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.42      0.27    10.91    11.96 1.00     3301     2477

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

fit1 plots

plot(fit1, N = 3, ask=FALSE)

fit 1 predictions

fit1.prediction <- cbind(datbcd, pred1=predict(fit1)[,"Estimate"])

fit1.prediction %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm), alpha=.1) +
  geom_line(aes(y=pred1), color="blue") 

fit 1 plot predictions per pop

  • Please write code to generate the plot below (you can do this by adding one function call to the code on the previous slide)
  • In your own words, how did it do?

Allow a different Hmax for each pop

  • The plants from different pops grow to different heights. Maybe we should model that.
  • With only 3 pops we will use a fixed effect. (Need 5+ for random effect)
  • See below for how we specify that Hmax should take on a different value for each pop
f2 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + k + delta ~ 1,
                  Hmax ~ pop,
                  nl=TRUE)

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b     popCC                  Hmax       (vectorized)
               (flat)     b    popDPR                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)

Set the priors

  • Do these seem reasonable to you?
prior2 <- c(set_prior("normal(5,3)", nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta")
)

Fit model 2

fit2 <- brm(formula=f2, data=datbcd, prior=prior2, cores = 4)
Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1

Fit 2 summary

Any Problems?

summary(fit2)
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         k ~ 1
         delta ~ 1
         Hmax ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      4.02      0.55     2.90     5.09 1.00     3269     2623
k_Intercept         0.66      0.01     0.64     0.68 1.00     2303     2101
delta_Intercept     6.04      0.45     5.20     6.95 1.00     2511     2354
Hmax_Intercept     65.52      1.53    62.70    68.74 1.00     2437     1920
Hmax_popCC         10.42     11.66   -12.15    33.50 1.00     3890     2734
Hmax_popDPR        48.86     12.18    25.47    72.77 1.00     3424     2261

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.28      0.26    10.79    11.80 1.00     4021     3007

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

Fit 2 plots

Any problems?

plot(fit2, N = 3, ask=FALSE)

Fit 2 plot the predictions

datbcd$pred2 <- predict(fit2)[,"Estimate"]

datbcd %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop), alpha=.3) +
  geom_line(aes(y=pred2, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2", guide="none") +
  facet_wrap(~pop)

Fit 2 Extend the CC and DPR curves out for all days

  • It would be nice to get predictions for CC and DPR all the way out to day 200
  • To do this, we create a data frame containing each day and pop for which we would like a prediction
  • The function expand_grid creates a data frame with all combinations of the items it is given (days and pops in this case).
pred.df <- expand_grid(day=min(datbcd$day):max(datbcd$day),
                       pop=unique(datbcd$pop))
pred.df
# A tibble: 537 × 2
     day pop  
   <int> <chr>
 1    30 CC   
 2    30 BH   
 3    30 DPR  
 4    31 CC   
 5    31 BH   
 6    31 DPR  
 7    32 CC   
 8    32 BH   
 9    32 DPR  
10    33 CC   
# ℹ 527 more rows

Fit 2 Extend the CC and DPR curves out for all days (slide 2)

  • Now we use our prediction data frame to generate predictions
fit2.predictions <- cbind(pred.df, prediction=predict(fit2, newdata = pred.df)[,"Estimate"]) %>%
  full_join(datbcd, by=c("day", "pop"))

fit2.predictions
     day pop prediction plantID survey_date block row  col   mf rep height_cm
1     30  CC   3.903193    D14B  2023-01-27    D1   4    B    5  12       3.7
2     30  CC   3.903193    D16B  2023-01-27    D1   6    B    3  10       2.3
3     30  CC   3.903193   D112A  2023-01-27    D1  12    A    7   4       4.1
4     30  CC   3.903193   D228B  2023-01-27    D2  28    B    2  13       4.1
5     30  CC   3.903193   D235A  2023-01-27    D2  35    A    1  11       5.4
6     30  CC   3.903193   D239A  2023-01-27    D2  39    A    5   7       3.6
7     30  CC   3.903193   D239B  2023-01-27    D2  39    B    4   5       5.1
8     30  CC   3.903193    D14C  2023-01-27    D1   4    C    4   9       0.9
9     30  CC   3.903193   F112A  2023-01-27    F1  12    A    6   3       6.5
10    30  CC   3.903193   F120B  2023-01-27    F1  20    B    2   4       6.0
11    30  CC   3.903193   F231A  2023-01-27    F2  31    A    7   8       4.9
12    30  CC   3.903193   F117D  2023-01-27    F1  17    D    3  13       4.7
13    30  CC   3.903193   F224C  2023-01-27    F2  24    C    1   6       3.4
14    30  CC   3.903193   F234D  2023-01-27    F2  34    D    3  12       5.6
15    30  CC   3.903193   H111B  2023-01-27    H1  11    B    7   2       6.8
16    30  CC   3.903193   H115B  2023-01-27    H1  15    B    3   3       5.5
17    30  CC   3.903193   H116A  2023-01-27    H1  16    A    4  10       6.3
18    30  CC   3.903193   H116B  2023-01-27    H1  16    B    6   6       4.4
19    30  CC   3.903193   H221A  2023-01-27    H2  21    A    3   5       4.8
20    30  CC   3.903193   H234B  2023-01-27    H2  34    B    4   7       3.4
21    30  CC   3.903193    H19D  2023-01-27    H1   9    D    2   2       4.8
22    30  CC   3.903193   H236C  2023-01-27    H2  36    C    2   1       5.6
23    30  CC   3.903193    J19D  2023-01-27    J1   9    D    2   6       3.8
24    30  CC   3.903193   J116C  2023-01-27    J1  16    C    6   4       5.6
25    30  CC   3.903193   J236C  2023-01-27    J2  36    C    5   5       4.2
26    30  CC   3.903193   J239C  2023-01-27    J2  39    C    2   3       6.9
27    30  CC   3.903193   J110A  2023-01-27    J1  10    A    3   2       6.3
28    30  CC   3.903193   J112B  2023-01-27    J1  12    B    4   4       4.1
29    30  CC   3.903193   J238B  2023-01-27    J2  38    B    3   9       4.6
30    30  CC   3.903193    L18A  2023-01-27    L1   8    A    4   6       3.4
31    30  CC   3.903193   L232A  2023-01-27    L2  32    A    5   3       5.5
32    30  CC   3.903193   L114D  2023-01-27    L1  14    D    1   1       3.9
33    30  CC   3.903193   L117C  2023-01-27    L1  17    C    7   3       4.2
34    30  BH   4.016931    D15B  2023-01-27    D1   5    B    5  24       3.9
35    30  BH   4.016931    D18B  2023-01-27    D1   8    B    6   7       2.2
36    30  BH   4.016931   D111B  2023-01-27    D1  11    B    7   4       3.7
37    30  BH   4.016931   D119B  2023-01-27    D1  19    B    2   8       3.7
38    30  BH   4.016931   D120A  2023-01-27    D1  20    A    4   2       3.3
39    30  BH   4.016931   D226B  2023-01-27    D2  26    B    5  15       2.7
40    30  BH   4.016931   D230B  2023-01-27    D2  30    B    4   3       3.0
41    30  BH   4.016931   D233B  2023-01-27    D2  33    B    6   8       2.2
42    30  BH   4.016931   D235B  2023-01-27    D2  35    B    5  25       2.6
43    30  BH   4.016931    D16C  2023-01-27    D1   6    C    4  12       3.1
44    30  BH   4.016931    D18D  2023-01-27    D1   8    D    7  17       1.6
45    30  BH   4.016931    D19C  2023-01-27    D1   9    C    1   8       2.1
46    30  BH   4.016931   D117C  2023-01-27    D1  17    C    5  14       3.2
47    30  BH   4.016931   D118D  2023-01-27    D1  18    D    3   5       2.3
48    30  BH   4.016931   D119D  2023-01-27    D1  19    D    5   9       1.4
49    30  BH   4.016931   D222C  2023-01-27    D2  22    C    7  18       1.4
50    30  BH   4.016931   D227C  2023-01-27    D2  27    C    7   5       3.1
51    30  BH   4.016931   D229C  2023-01-27    D2  29    C    4  13       3.3
52    30  BH   4.016931   D229D  2023-01-27    D2  29    D    2   9       2.9
53    30  BH   4.016931   D230D  2023-01-27    D2  30    D    1  10       2.3
54    30  BH   4.016931   D235D  2023-01-27    D2  35    D    3   6       3.7
55    30  BH   4.016931   F110A  2023-01-27    F1  10    A    7  19       2.9
56    30  BH   4.016931   F115A  2023-01-27    F1  15    A    6   9       3.6
57    30  BH   4.016931   F117B  2023-01-27    F1  17    B    2  10       3.4
58    30  BH   4.016931   F227A  2023-01-27    F2  27    A    7   7       3.4
59    30  BH   4.016931   F229B  2023-01-27    F2  29    B    4   5       4.8
60    30  BH   4.016931   F240A  2023-01-27    F2  40    A    3   8       3.6
61    30  BH   4.016931   F240B  2023-01-27    F2  40    B    5  12       6.5
62    30  BH   4.016931    F13C  2023-01-27    F1   3    C    5  16       2.9
63    30  BH   4.016931    F15C  2023-01-27    F1   5    C    5  26       2.9
64    30  BH   4.016931    F17D  2023-01-27    F1   7    D    7   6       3.5
65    30  BH   4.016931   F113C  2023-01-27    F1  13    C    3   7       2.8
66    30  BH   4.016931   F115C  2023-01-27    F1  15    C    4   4       3.4
67    30  BH   4.016931   F120C  2023-01-27    F1  20    C    5  11       5.1
68    30  BH   4.016931   F120D  2023-01-27    F1  20    D    5   1       2.9
69    30  BH   4.016931   F225C  2023-01-27    F2  25    C    7  20       2.4
70    30  BH   4.016931   F231C  2023-01-27    F2  31    C    6  10       3.0
71    30  BH   4.016931   F233C  2023-01-27    F2  33    C    5  17       3.0
72    30  BH   4.016931   F233D  2023-01-27    F2  33    D    2  11       3.6
73    30  BH   4.016931   F235C  2023-01-27    F2  35    C    5  27       3.2
74    30  BH   4.016931   F235D  2023-01-27    F2  35    D    2   1       3.7
75    30  BH   4.016931    H14B  2023-01-27    H1   4    B    2  12       4.5
76    30  BH   4.016931    H18B  2023-01-27    H1   8    B    3   9       3.2
77    30  BH   4.016931   H114A  2023-01-27    H1  14    A    5   3       6.9
78    30  BH   4.016931   H117A  2023-01-27    H1  17    A    4   6       2.3
79    30  BH   4.016931   H118B  2023-01-27    H1  18    B    5  18       3.5
80    30  BH   4.016931   H224B  2023-01-27    H2  24    B    5  29       4.5
81    30  BH   4.016931   H227A  2023-01-27    H2  27    A    5   4       4.5
82    30  BH   4.016931   H229B  2023-01-27    H2  29    B    5  19       3.2
83    30  BH   4.016931   H233A  2023-01-27    H2  33    A    3  10       2.5
84    30  BH   4.016931   H233B  2023-01-27    H2  33    B    4   7       2.9
85    30  BH   4.016931   H239A  2023-01-27    H2  39    A    7   9       3.0
86    30  BH   4.016931   H239B  2023-01-27    H2  39    B    6   2       3.5
87    30  BH   4.016931    H15C  2023-01-27    H1   5    C    5  28       4.6
88    30  BH   4.016931   H110C  2023-01-27    H1  10    C    6  11       4.1
89    30  BH   4.016931   H115D  2023-01-27    H1  15    D    5  13       4.1
90    30  BH   4.016931   H117C  2023-01-27    H1  17    C    2   2       3.5
91    30  BH   4.016931   H119C  2023-01-27    H1  19    C    7  21       4.5
92    30  BH   4.016931   H222D  2023-01-27    H2  22    D    2   3       2.8
93    30  BH   4.016931   H223C  2023-01-27    H2  23    C    7  22       2.8
94    30  BH   4.016931   H227D  2023-01-27    H2  27    D    6  12       2.8
95    30  BH   4.016931   H235C  2023-01-27    H2  35    C    2  13       2.5
96    30  BH   4.016931    J19C  2023-01-27    J1   9    C    5  30       4.4
97    30  BH   4.016931   J110C  2023-01-27    J1  10    C    3   1       1.7
98    30  BH   4.016931   J112C  2023-01-27    J1  12    C    7  23       2.0
99    30  BH   4.016931   J114D  2023-01-27    J1  14    D    4   8       3.3
100   30  BH   4.016931   J116D  2023-01-27    J1  16    D    5   5       2.7
101   30  BH   4.016931   J119C  2023-01-27    J1  19    C    5  20       4.3
102   30  BH   4.016931   J223D  2023-01-27    J2  23    D    4   9       4.3
103   30  BH   4.016931   J225D  2023-01-27    J2  25    D    2   5       1.9
104   30  BH   4.016931   J230C  2023-01-27    J2  30    C    4  14       4.4
105   30  BH   4.016931   J231C  2023-01-27    J2  31    C    5   6       6.1
106   30  BH   4.016931   J231D  2023-01-27    J2  31    D    7  14       2.7
107   30  BH   4.016931   J232D  2023-01-27    J2  32    D    3  12       3.5
108   30  BH   4.016931   J235C  2023-01-27    J2  35    C    7  11       1.5
109   30  BH   4.016931   J236D  2023-01-27    J2  36    D    5  21       3.0
110   30  BH   4.016931    J14B  2023-01-27    J1   4    B    2   4       2.3
111   30  BH   4.016931   J119B  2023-01-27    J1  19    B    3  11       3.8
112   30  BH   4.016931   J120B  2023-01-27    J1  20    B    7  10       2.2
113   30  BH   4.016931   J241B  2023-01-27    J2  41    B    6   4       5.5
114   30  BH   4.016931   J233A  2023-01-27    J2  33    A    3   2       3.0
115   30  BH   4.016931   J237B  2023-01-27    J2  37    B    1   7       3.0
116   30  BH   4.016931    L15B  2023-01-27    L1   5    B    2   6       2.8
117   30  BH   4.016931    L18B  2023-01-27    L1   8    B    7   2       3.3
118   30  BH   4.016931   L110A  2023-01-27    L1  10    A    2  14       1.4
119   30  BH   4.016931   L111B  2023-01-27    L1  11    B    3   3       2.7
120   30  BH   4.016931   L115A  2023-01-27    L1  15    A    5   7       3.0
121   30  BH   4.016931   L120A  2023-01-27    L1  20    A    7  15       2.3
122   30  BH   4.016931    L14D  2023-01-27    L1   4    D    4  10       2.1
123   30  BH   4.016931   L113C  2023-01-27    L1  13    C    3  13       2.4
124   30  BH   4.016931   L113D  2023-01-27    L1  13    D    5  22       3.7
125   30  BH   4.016931   L116C  2023-01-27    L1  16    C    6   5       4.4
126   30  BH   4.016931   L118D  2023-01-27    L1  18    D    7  12       3.7
127   30  BH   4.016931   L228C  2023-01-27    L2  28    C    6   6       3.1
128   30  BH   4.016931   L228D  2023-01-27    L2  28    D    2   7       2.9
129   30  BH   4.016931   L229D  2023-01-27    L2  29    D    3  14       1.2
130   30  BH   4.016931   L230C  2023-01-27    L2  30    C    4  11       4.5
131   30  BH   4.016931   L234C  2023-01-27    L2  34    C    5   8       3.6
132   30  BH   4.016931   L237D  2023-01-27    L2  37    D    7  16       2.6
133   30  BH   4.016931   L240D  2023-01-27    L2  40    D    7  13       3.8
134   30  BH   4.016931   L241D  2023-01-27    L2  41    D    4   1       3.3
135   30  BH   4.016931   L242C  2023-01-27    L2  42    C    7   3       2.0
136   30 DPR   4.109382   D117A  2023-01-27    D1  17    A    6   4       4.1
137   30 DPR   4.109382   D225A  2023-01-27    D2  25    A    2   7       8.5
138   30 DPR   4.109382   D227B  2023-01-27    D2  27    B    4   9       7.0
139   30 DPR   4.109382    F19A  2023-01-27    F1   9    A    2   5       6.5
140   30 DPR   4.109382   F232A  2023-01-27    F2  32    A    6   2       2.5
141   30 DPR   4.109382   F241A  2023-01-27    F2  41    A    2   4       5.5
142   30 DPR   4.109382    F16C  2023-01-27    F1   6    C    6   9       2.5
143   30 DPR   4.109382   F240C  2023-01-27    F2  40    C    7   4       5.4
144   30 DPR   4.109382   H117B  2023-01-27    H1  17    B    5   2       6.4
145   30 DPR   4.109382   H226A  2023-01-27    H2  26    A    4  13       3.5
146   30 DPR   4.109382    H17C  2023-01-27    H1   7    C    2  10       7.8
147   30 DPR   4.109382   H225C  2023-01-27    H2  25    C    2   6       7.5
148   30 DPR   4.109382   J121D  2023-01-27    J1  21    D    4   8       3.9
149   30 DPR   4.109382    J16B  2023-01-27    J1   6    B    6  10       4.1
150   30 DPR   4.109382   J232A  2023-01-27    J2  32    A    2  12       7.5
151   30 DPR   4.109382   J232B  2023-01-27    J2  32    B    6  11       3.2
152   30 DPR   4.109382   L119B  2023-01-27    L1  19    B    4  11       4.3
153   30 DPR   4.109382   L237B  2023-01-27    L2  37    B    6   3       2.3
154   30 DPR   4.109382   L116D  2023-01-27    L1  16    D    6  13       4.5
155   30 DPR   4.109382   L233C  2023-01-27    L2  33    C    2   2       5.4
156   31  CC   4.137701    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
157   31  BH   4.165617    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
158   31 DPR   4.131556    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
159   32  CC   4.360065    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
160   32  BH   4.006778    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
161   32 DPR   4.146713    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
162   33  CC   4.208230    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
163   33  BH   3.828360    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
164   33 DPR   4.221002    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
165   34  CC   4.110044    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
166   34  BH   3.984476    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
167   34 DPR   3.781387    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
168   35  CC   4.242742    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
169   35  BH   4.298571    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
170   35 DPR   3.768074    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
171   36  CC   4.146959    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
172   36  BH   4.008429    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
173   36 DPR   3.966143    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
174   37  CC   4.041823    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
175   37  BH   3.853087    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
176   37 DPR   3.961463    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
177   38  CC   3.992897    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
178   38  BH   4.036273    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
179   38 DPR   4.188412    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
180   39  CC   4.427747    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
181   39  BH   4.127631    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
182   39 DPR   3.880306    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
183   40  CC   3.586425    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
184   40  BH   4.150061    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
185   40 DPR   4.029894    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
186   41  CC   4.205030    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
187   41  BH   4.080003    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
188   41 DPR   4.214745    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
189   42  CC   4.323083    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
190   42  BH   3.964821    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
191   42 DPR   4.107392    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
192   43  CC   4.050626    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
193   43  BH   4.106451    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
194   43 DPR   4.062066    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
195   44  CC   4.020870    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
196   44  BH   4.026386    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
197   44 DPR   4.077503    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
198   45  CC   3.723449    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
199   45  BH   4.321398    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
200   45 DPR   3.977987    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
201   46  CC   4.296484    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
202   46  BH   3.734842    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
203   46 DPR   3.819111    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
204   47  CC   3.888357    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
205   47  BH   4.134525    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
206   47 DPR   4.206847    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
207   48  CC   4.103067    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
208   48  BH   4.148949    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
209   48 DPR   4.103073    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
210   49  CC   3.861544    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
211   49  BH   4.096552    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
212   49 DPR   4.222690    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
213   50  CC   4.240629    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
214   50  BH   4.289048    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
215   50 DPR   4.112767    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
216   51  CC   4.221596    D14B  2023-02-17    D1   4    B    5  12       3.4
217   51  CC   4.221596   D112A  2023-02-17    D1  12    A    7   4       4.3
218   51  CC   4.221596   D228B  2023-02-17    D2  28    B    2  13       4.0
219   51  CC   4.221596   D235A  2023-02-17    D2  35    A    1  11       5.5
220   51  CC   4.221596   D239A  2023-02-17    D2  39    A    5   7       3.8
221   51  CC   4.221596   D239B  2023-02-17    D2  39    B    4   5       5.2
222   51  CC   4.221596   F112A  2023-02-17    F1  12    A    6   3       6.1
223   51  CC   4.221596   F120B  2023-02-17    F1  20    B    2   4       5.6
224   51  CC   4.221596   F231A  2023-02-17    F2  31    A    7   8       4.1
225   51  CC   4.221596   F117D  2023-02-17    F1  17    D    3  13       4.6
226   51  CC   4.221596   F234D  2023-02-17    F2  34    D    3  12       5.5
227   51  CC   4.221596   H111B  2023-02-17    H1  11    B    7   2       7.6
228   51  CC   4.221596   H115B  2023-02-17    H1  15    B    3   3       6.1
229   51  CC   4.221596   H116A  2023-02-17    H1  16    A    4  10       6.1
230   51  CC   4.221596   H234B  2023-02-17    H2  34    B    4   7       2.9
231   51  CC   4.221596    H19D  2023-02-17    H1   9    D    2   2       4.0
232   51  CC   4.221596   H236C  2023-02-17    H2  36    C    2   1       6.2
233   51  CC   4.221596   J110A  2023-02-17    J1  10    A    3   2       6.5
234   51  CC   4.221596   J112B  2023-02-17    J1  12    B    4   4       4.5
235   51  CC   4.221596   J228A  2023-02-17    J2  28    A    5   6       3.6
236   51  CC   4.221596   J238B  2023-02-17    J2  38    B    3   9       3.5
237   51  CC   4.221596    J19D  2023-02-17    J1   9    D    2   6       3.1
238   51  CC   4.221596   J236C  2023-02-17    J2  36    C    5   5       4.9
239   51  CC   4.221596   J239C  2023-02-17    J2  39    C    2   3       7.7
240   51  CC   4.221596    L18A  2023-02-17    L1   8    A    4   6       4.1
241   51  CC   4.221596   L232A  2023-02-17    L2  32    A    5   3       5.9
242   51  CC   4.221596   L117C  2023-02-17    L1  17    C    7   3       4.9
243   51  BH   4.066216    D15B  2023-02-17    D1   5    B    5  24       4.3
244   51  BH   4.066216    D18B  2023-02-17    D1   8    B    6   7       2.9
245   51  BH   4.066216   D111B  2023-02-17    D1  11    B    7   4       3.5
246   51  BH   4.066216   D119B  2023-02-17    D1  19    B    2   8       3.3
247   51  BH   4.066216   D120A  2023-02-17    D1  20    A    4   2       3.9
248   51  BH   4.066216   D226B  2023-02-17    D2  26    B    5  15       3.0
249   51  BH   4.066216   D230B  2023-02-17    D2  30    B    4   3       3.9
250   51  BH   4.066216   D233B  2023-02-17    D2  33    B    6   8       2.2
251   51  BH   4.066216   D235B  2023-02-17    D2  35    B    5  25       1.9
252   51  BH   4.066216    D16C  2023-02-17    D1   6    C    4  12       4.7
253   51  BH   4.066216    D18D  2023-02-17    D1   8    D    7  17       1.5
254   51  BH   4.066216    D19C  2023-02-17    D1   9    C    1   8      24.0
255   51  BH   4.066216   D117C  2023-02-17    D1  17    C    5  14       2.9
256   51  BH   4.066216   D118D  2023-02-17    D1  18    D    3   5       2.2
257   51  BH   4.066216   D119D  2023-02-17    D1  19    D    5   9       0.9
258   51  BH   4.066216   D227C  2023-02-17    D2  27    C    7   5       3.4
259   51  BH   4.066216   D229C  2023-02-17    D2  29    C    4  13       3.4
260   51  BH   4.066216   D229D  2023-02-17    D2  29    D    2   9       3.3
261   51  BH   4.066216   D230D  2023-02-17    D2  30    D    1  10       2.1
262   51  BH   4.066216   D235D  2023-02-17    D2  35    D    3   6       4.0
263   51  BH   4.066216   F115A  2023-02-17    F1  15    A    6   9       4.1
264   51  BH   4.066216   F117B  2023-02-17    F1  17    B    2  10       3.0
265   51  BH   4.066216   F229B  2023-02-17    F2  29    B    4   5       5.1
266   51  BH   4.066216   F240B  2023-02-17    F2  40    B    5  12       4.2
267   51  BH   4.066216    F15C  2023-02-17    F1   5    C    5  26       2.9
268   51  BH   4.066216    F17D  2023-02-17    F1   7    D    7   6       3.1
269   51  BH   4.066216   F113C  2023-02-17    F1  13    C    3   7       2.9
270   51  BH   4.066216   F115C  2023-02-17    F1  15    C    4   4       3.5
271   51  BH   4.066216   F120D  2023-02-17    F1  20    D    5   1       1.7
272   51  BH   4.066216   F225C  2023-02-17    F2  25    C    7  20       2.2
273   51  BH   4.066216   F233C  2023-02-17    F2  33    C    5  17       3.1
274   51  BH   4.066216   F233D  2023-02-17    F2  33    D    2  11       3.7
275   51  BH   4.066216   F235C  2023-02-17    F2  35    C    5  27       3.5
276   51  BH   4.066216   F235D  2023-02-17    F2  35    D    2   1       3.6
277   51  BH   4.066216    H14B  2023-02-17    H1   4    B    2  12       5.2
278   51  BH   4.066216    H18B  2023-02-17    H1   8    B    3   9       3.3
279   51  BH   4.066216   H114A  2023-02-17    H1  14    A    5   3       7.4
280   51  BH   4.066216   H117A  2023-02-17    H1  17    A    4   6       3.5
281   51  BH   4.066216   H118B  2023-02-17    H1  18    B    5  18       3.2
282   51  BH   4.066216   H224B  2023-02-17    H2  24    B    5  29       4.6
283   51  BH   4.066216   H227A  2023-02-17    H2  27    A    5   4       3.9
284   51  BH   4.066216   H229B  2023-02-17    H2  29    B    5  19       3.1
285   51  BH   4.066216   H233B  2023-02-17    H2  33    B    4   7       2.5
286   51  BH   4.066216   H239B  2023-02-17    H2  39    B    6   2       3.9
287   51  BH   4.066216    H15C  2023-02-17    H1   5    C    5  28       4.7
288   51  BH   4.066216   H110C  2023-02-17    H1  10    C    6  11       4.0
289   51  BH   4.066216   H115D  2023-02-17    H1  15    D    5  13       4.5
290   51  BH   4.066216   H119C  2023-02-17    H1  19    C    7  21       4.1
291   51  BH   4.066216   H222D  2023-02-17    H2  22    D    2   3       2.4
292   51  BH   4.066216   H223C  2023-02-17    H2  23    C    7  22       3.2
293   51  BH   4.066216   H227D  2023-02-17    H2  27    D    6  12       3.8
294   51  BH   4.066216   H235C  2023-02-17    H2  35    C    2  13       2.2
295   51  BH   4.066216    J14B  2023-02-17    J1   4    B    2   4       2.6
296   51  BH   4.066216   J119B  2023-02-17    J1  19    B    3  11       4.2
297   51  BH   4.066216   J233A  2023-02-17    J2  33    A    3   2       3.1
298   51  BH   4.066216   J237B  2023-02-17    J2  37    B    1   7       2.7
299   51  BH   4.066216   J241B  2023-02-17    J2  41    B    6   4       4.5
300   51  BH   4.066216    J19C  2023-02-17    J1   9    C    5  30       5.5
301   51  BH   4.066216   J110C  2023-02-17    J1  10    C    3   1       1.9
302   51  BH   4.066216   J114D  2023-02-17    J1  14    D    4   8       3.4
303   51  BH   4.066216   J119C  2023-02-17    J1  19    C    5  20       5.3
304   51  BH   4.066216   J223D  2023-02-17    J2  23    D    4   9       1.7
305   51  BH   4.066216   J225D  2023-02-17    J2  25    D    2   5       1.9
306   51  BH   4.066216   J230C  2023-02-17    J2  30    C    4  14       5.3
307   51  BH   4.066216   J231C  2023-02-17    J2  31    C    5   6       7.1
308   51  BH   4.066216   J231D  2023-02-17    J2  31    D    7  14       4.8
309   51  BH   4.066216   J232D  2023-02-17    J2  32    D    3  12       3.4
310   51  BH   4.066216   J236D  2023-02-17    J2  36    D    5  21      25.0
311   51  BH   4.066216    L15B  2023-02-17    L1   5    B    2   6       2.9
312   51  BH   4.066216    L18B  2023-02-17    L1   8    B    7   2       2.7
313   51  BH   4.066216   L110A  2023-02-17    L1  10    A    2  14       0.6
314   51  BH   4.066216   L111B  2023-02-17    L1  11    B    3   3       2.1
315   51  BH   4.066216   L120A  2023-02-17    L1  20    A    7  15       2.3
316   51  BH   4.066216    L14D  2023-02-17    L1   4    D    4  10       2.2
317   51  BH   4.066216   L113C  2023-02-17    L1  13    C    3  13       2.2
318   51  BH   4.066216   L113D  2023-02-17    L1  13    D    5  22       4.7
319   51  BH   4.066216   L116C  2023-02-17    L1  16    C    6   5       4.4
320   51  BH   4.066216   L228C  2023-02-17    L2  28    C    6   6       2.4
321   51  BH   4.066216   L228D  2023-02-17    L2  28    D    2   7       2.8
322   51  BH   4.066216   L229D  2023-02-17    L2  29    D    3  14       1.9
323   51  BH   4.066216   L230C  2023-02-17    L2  30    C    4  11       4.4
324   51  BH   4.066216   L234C  2023-02-17    L2  34    C    5   8       4.1
325   51  BH   4.066216   L237D  2023-02-17    L2  37    D    7  16       3.6
326   51  BH   4.066216   L240D  2023-02-17    L2  40    D    7  13       4.1
327   51  BH   4.066216   L241D  2023-02-17    L2  41    D    4   1       3.5
328   51 DPR   4.009491   D225A  2023-02-17    D2  25    A    2   7       8.9
329   51 DPR   4.009491   D227B  2023-02-17    D2  27    B    4   9       7.2
330   51 DPR   4.009491    F19A  2023-02-17    F1   9    A    2   5       6.0
331   51 DPR   4.009491   F232A  2023-02-17    F2  32    A    6   2       2.0
332   51 DPR   4.009491    F16C  2023-02-17    F1   6    C    6   9       2.6
333   51 DPR   4.009491   F240C  2023-02-17    F2  40    C    7   4       6.3
334   51 DPR   4.009491   H117B  2023-02-17    H1  17    B    5   2       8.1
335   51 DPR   4.009491    H17C  2023-02-17    H1   7    C    2  10       8.2
336   51 DPR   4.009491   H225C  2023-02-17    H2  25    C    2   6       7.8
337   51 DPR   4.009491    J16B  2023-02-17    J1   6    B    6  10       3.6
338   51 DPR   4.009491   J232B  2023-02-17    J2  32    B    6  11       2.1
339   51 DPR   4.009491   L119B  2023-02-17    L1  19    B    4  11       3.7
340   51 DPR   4.009491   L237B  2023-02-17    L2  37    B    6   3       4.0
341   51 DPR   4.009491   L116D  2023-02-17    L1  16    D    6  13       4.0
342   51 DPR   4.009491   L233C  2023-02-17    L2  33    C    2   2       5.7
343   52  CC   4.240536    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
344   52  BH   4.312868    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
345   52 DPR   3.971057    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
346   53  CC   4.034819    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
347   53  BH   4.388159    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
348   53 DPR   4.268878    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
349   54  CC   4.144131    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
350   54  BH   4.563911    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
351   54 DPR   4.644620    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
352   55  CC   4.207493    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
353   55  BH   4.190144    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
354   55 DPR   4.573265    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
355   56  CC   4.251412    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
356   56  BH   4.255124    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
357   56 DPR   4.180665    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
358   57  CC   4.453380    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
359   57  BH   4.361777    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
360   57 DPR   4.641373    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
361   58  CC   4.357757    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
362   58  BH   3.966031    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
363   58 DPR   4.297890    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
364   59  CC   4.293584    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
365   59  BH   4.494128    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
366   59 DPR   4.266290    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
367   60  CC   4.268374    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
368   60  BH   4.435386    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
369   60 DPR   4.568954    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
370   61  CC   4.279176    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
371   61  BH   4.433645    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
372   61 DPR   4.390738    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
373   62  CC   4.360378    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
374   62  BH   4.290543    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
375   62 DPR   4.365633    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
376   63  CC   4.578521    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
377   63  BH   4.317363    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
378   63 DPR   5.013411    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
379   64  CC   4.469379    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
380   64  BH   4.203130    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
381   64 DPR   4.434334    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
382   65  CC   4.275072   D112A  2023-03-03    D1  12    A    7   4       5.7
383   65  CC   4.275072   D235A  2023-03-03    D2  35    A    1  11       5.9
384   65  CC   4.275072   D239A  2023-03-03    D2  39    A    5   7       3.9
385   65  CC   4.275072   D239B  2023-03-03    D2  39    B    4   5       5.8
386   65  CC   4.275072   F231A  2023-03-03    F2  31    A    7   8       5.4
387   65  CC   4.275072   F117D  2023-03-03    F1  17    D    3  13       5.7
388   65  CC   4.275072   H111B  2023-03-03    H1  11    B    7   2       8.6
389   65  CC   4.275072   H115B  2023-03-03    H1  15    B    3   3       6.4
390   65  CC   4.275072   H116A  2023-03-03    H1  16    A    4  10       8.0
391   65  CC   4.275072   H234B  2023-03-03    H2  34    B    4   7       3.5
392   65  CC   4.275072   H236C  2023-03-03    H2  36    C    2   1       5.9
393   65  CC   4.275072   J112B  2023-03-03    J1  12    B    4   4       5.8
394   65  CC   4.275072   J228A  2023-03-03    J2  28    A    5   6       4.3
395   65  CC   4.275072   J236C  2023-03-03    J2  36    C    5   5       5.2
396   65  CC   4.275072   J239C  2023-03-03    J2  39    C    2   3       8.2
397   65  CC   4.275072    L18A  2023-03-03    L1   8    A    4   6       5.5
398   65  CC   4.275072   L232A  2023-03-03    L2  32    A    5   3       6.3
399   65  CC   4.275072   L117C  2023-03-03    L1  17    C    7   3       4.8
400   65  BH   4.775040    D18B  2023-03-03    D1   8    B    6   7       3.3
401   65  BH   4.775040   D120A  2023-03-03    D1  20    A    4   2       6.5
402   65  BH   4.775040   D226B  2023-03-03    D2  26    B    5  15       4.2
403   65  BH   4.775040   D230B  2023-03-03    D2  30    B    4   3       4.9
404   65  BH   4.775040   D233B  2023-03-03    D2  33    B    6   8       4.0
405   65  BH   4.775040    D16C  2023-03-03    D1   6    C    4  12       5.9
406   65  BH   4.775040    D19C  2023-03-03    D1   9    C    1   8       2.7
407   65  BH   4.775040   D119D  2023-03-03    D1  19    D    5   9       1.4
408   65  BH   4.775040   D229C  2023-03-03    D2  29    C    4  13       4.7
409   65  BH   4.775040   D229D  2023-03-03    D2  29    D    2   9       4.8
410   65  BH   4.775040   D235D  2023-03-03    D2  35    D    3   6       5.2
411   65  BH   4.775040   F115A  2023-03-03    F1  15    A    6   9       5.5
412   65  BH   4.775040   F229B  2023-03-03    F2  29    B    4   5       5.0
413   65  BH   4.775040   F240B  2023-03-03    F2  40    B    5  12       4.3
414   65  BH   4.775040    F15C  2023-03-03    F1   5    C    5  26       2.9
415   65  BH   4.775040   F113C  2023-03-03    F1  13    C    3   7       3.8
416   65  BH   4.775040   F115C  2023-03-03    F1  15    C    4   4       4.9
417   65  BH   4.775040   F235C  2023-03-03    F2  35    C    5  27       4.2
418   65  BH   4.775040   F235D  2023-03-03    F2  35    D    2   1       4.5
419   65  BH   4.775040    H14B  2023-03-03    H1   4    B    2  12       5.1
420   65  BH   4.775040    H18B  2023-03-03    H1   8    B    3   9       3.5
421   65  BH   4.775040   H114A  2023-03-03    H1  14    A    5   3       6.9
422   65  BH   4.775040   H117A  2023-03-03    H1  17    A    4   6       3.3
423   65  BH   4.775040   H118B  2023-03-03    H1  18    B    5  18       5.0
424   65  BH   4.775040   H224B  2023-03-03    H2  24    B    5  29       9.3
425   65  BH   4.775040   H227A  2023-03-03    H2  27    A    5   4       3.9
426   65  BH   4.775040   H229B  2023-03-03    H2  29    B    5  19       4.4
427   65  BH   4.775040   H233B  2023-03-03    H2  33    B    4   7       2.5
428   65  BH   4.775040   H239B  2023-03-03    H2  39    B    6   2       3.5
429   65  BH   4.775040    H15C  2023-03-03    H1   5    C    5  28       5.3
430   65  BH   4.775040   H110C  2023-03-03    H1  10    C    6  11       3.7
431   65  BH   4.775040   H115D  2023-03-03    H1  15    D    5  13       4.8
432   65  BH   4.775040   H119C  2023-03-03    H1  19    C    7  21       4.6
433   65  BH   4.775040   H222D  2023-03-03    H2  22    D    2   3       3.1
434   65  BH   4.775040   H223C  2023-03-03    H2  23    C    7  22       3.4
435   65  BH   4.775040   H235C  2023-03-03    H2  35    C    2  13       2.8
436   65  BH   4.775040    J14B  2023-03-03    J1   4    B    2   4       3.6
437   65  BH   4.775040   J119B  2023-03-03    J1  19    B    3  11       5.6
438   65  BH   4.775040   J233A  2023-03-03    J2  33    A    3   2       4.6
439   65  BH   4.775040   J237B  2023-03-03    J2  37    B    1   7       4.8
440   65  BH   4.775040    J19C  2023-03-03    J1   9    C    5  30       5.6
441   65  BH   4.775040   J110C  2023-03-03    J1  10    C    3   1       2.2
442   65  BH   4.775040   J114D  2023-03-03    J1  14    D    4   8       3.5
443   65  BH   4.775040   J119C  2023-03-03    J1  19    C    5  20       5.2
444   65  BH   4.775040   J223D  2023-03-03    J2  23    D    4   9       4.9
445   65  BH   4.775040   J225D  2023-03-03    J2  25    D    2   5       1.9
446   65  BH   4.775040   J230C  2023-03-03    J2  30    C    4  14       6.3
447   65  BH   4.775040   J231C  2023-03-03    J2  31    C    5   6       7.8
448   65  BH   4.775040   J231D  2023-03-03    J2  31    D    7  14       7.7
449   65  BH   4.775040   J232D  2023-03-03    J2  32    D    3  12       3.5
450   65  BH   4.775040    L15B  2023-03-03    L1   5    B    2   6       3.8
451   65  BH   4.775040    L18B  2023-03-03    L1   8    B    7   2       2.3
452   65  BH   4.775040   L110A  2023-03-03    L1  10    A    2  14       1.2
453   65  BH   4.775040   L111B  2023-03-03    L1  11    B    3   3       2.5
454   65  BH   4.775040    L14D  2023-03-03    L1   4    D    4  10       3.6
455   65  BH   4.775040   L113C  2023-03-03    L1  13    C    3  13       2.4
456   65  BH   4.775040   L113D  2023-03-03    L1  13    D    5  22       4.6
457   65  BH   4.775040   L116C  2023-03-03    L1  16    C    6   5       4.8
458   65  BH   4.775040   L228C  2023-03-03    L2  28    C    6   6       2.2
459   65  BH   4.775040   L228D  2023-03-03    L2  28    D    2   7       3.5
460   65  BH   4.775040   L229D  2023-03-03    L2  29    D    3  14       2.1
461   65  BH   4.775040   L230C  2023-03-03    L2  30    C    4  11       5.1
462   65  BH   4.775040   L234C  2023-03-03    L2  34    C    5   8       4.4
463   65  BH   4.775040   L240D  2023-03-03    L2  40    D    7  13       4.3
464   65  BH   4.775040   L241D  2023-03-03    L2  41    D    4   1       3.2
465   65 DPR   4.577711   D225A  2023-03-03    D2  25    A    2   7       9.3
466   65 DPR   4.577711   D227B  2023-03-03    D2  27    B    4   9       7.9
467   65 DPR   4.577711   F232A  2023-03-03    F2  32    A    6   2       2.7
468   65 DPR   4.577711    F16C  2023-03-03    F1   6    C    6   9       4.3
469   65 DPR   4.577711   F240C  2023-03-03    F2  40    C    7   4       5.8
470   65 DPR   4.577711   H117B  2023-03-03    H1  17    B    5   2      10.1
471   65 DPR   4.577711    H17C  2023-03-03    H1   7    C    2  10       9.2
472   65 DPR   4.577711   H225C  2023-03-03    H2  25    C    2   6       7.6
473   65 DPR   4.577711    J16B  2023-03-03    J1   6    B    6  10       4.5
474   65 DPR   4.577711   J232B  2023-03-03    J2  32    B    6  11       4.4
475   65 DPR   4.577711   L237B  2023-03-03    L2  37    B    6   3       5.1
476   66  CC   4.462263    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
477   66  BH   4.707518    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
478   66 DPR   4.665711    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
479   67  CC   4.578495    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
480   67  BH   4.617266    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
481   67 DPR   4.941633    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
482   68  CC   4.773174    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
483   68  BH   5.088134    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
484   68 DPR   5.283460    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
485   69  CC   4.815656    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
486   69  BH   4.650572    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
487   69 DPR   5.149403    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
488   70  CC   4.647518    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
489   70  BH   4.593957    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
490   70 DPR   4.786918    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
491   71  CC   4.809500    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
492   71  BH   4.548636    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
493   71 DPR   5.410964    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
494   72  CC   4.960768    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
495   72  BH   4.551433    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
496   72 DPR   5.113534    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
497   73  CC   4.784553    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
498   73  BH   5.112026    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
499   73 DPR   5.480229    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
500   74  CC   4.967260    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
501   74  BH   4.890314    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
502   74 DPR   5.577908    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
503   75  CC   5.510598    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
504   75  BH   5.213073    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
505   75 DPR   5.535766    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
506   76  CC   5.332289    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
507   76  BH   5.101820    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
508   76 DPR   5.584991    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
509   77  CC   5.384407    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
510   77  BH   5.011949    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
511   77 DPR   5.868643    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
512   78  CC   5.791378    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
513   78  BH   5.131228    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
514   78 DPR   6.168828    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
515   79  CC   5.312003   D112A  2023-03-17    D1  12    A    7   4       5.4
516   79  CC   5.312003   D239A  2023-03-17    D2  39    A    5   7       3.9
517   79  CC   5.312003   D239B  2023-03-17    D2  39    B    4   5       5.2
518   79  CC   5.312003   F117D  2023-03-17    F1  17    D    3  13       7.2
519   79  CC   5.312003   H111B  2023-03-17    H1  11    B    7   2      12.3
520   79  CC   5.312003   H115B  2023-03-17    H1  15    B    3   3       9.7
521   79  CC   5.312003   H116A  2023-03-17    H1  16    A    4  10       9.1
522   79  CC   5.312003   H234B  2023-03-17    H2  34    B    4   7       3.4
523   79  CC   5.312003   H236C  2023-03-17    H2  36    C    2   1       6.1
524   79  CC   5.312003   J112B  2023-03-17    J1  12    B    4   4       6.3
525   79  CC   5.312003   J236C  2023-03-17    J2  36    C    5   5       7.2
526   79  CC   5.312003    L18A  2023-03-17    L1   8    A    4   6       4.2
527   79  CC   5.312003   L232A  2023-03-17    L2  32    A    5   3       5.9
528   79  CC   5.312003   L117C  2023-03-17    L1  17    C    7   3       5.4
529   79  BH   5.330013   D226B  2023-03-17    D2  26    B    5  15       5.3
530   79  BH   5.330013   D230B  2023-03-17    D2  30    B    4   3       6.1
531   79  BH   5.330013   D233B  2023-03-17    D2  33    B    6   8       5.6
532   79  BH   5.330013    D19C  2023-03-17    D1   9    C    1   8       2.6
533   79  BH   5.330013   D119D  2023-03-17    D1  19    D    5   9       2.1
534   79  BH   5.330013   D229C  2023-03-17    D2  29    C    4  13       4.3
535   79  BH   5.330013   D229D  2023-03-17    D2  29    D    2   9       4.9
536   79  BH   5.330013   D235D  2023-03-17    D2  35    D    3   6       5.3
537   79  BH   5.330013   F115A  2023-03-17    F1  15    A    6   9       8.0
538   79  BH   5.330013   F229B  2023-03-17    F2  29    B    4   5       5.9
539   79  BH   5.330013   F240B  2023-03-17    F2  40    B    5  12       5.1
540   79  BH   5.330013    F15C  2023-03-17    F1   5    C    5  26       3.0
541   79  BH   5.330013   F113C  2023-03-17    F1  13    C    3   7       4.1
542   79  BH   5.330013   F115C  2023-03-17    F1  15    C    4   4       4.5
543   79  BH   5.330013   F235C  2023-03-17    F2  35    C    5  27       5.5
544   79  BH   5.330013   F235D  2023-03-17    F2  35    D    2   1       4.6
545   79  BH   5.330013    H14B  2023-03-17    H1   4    B    2  12       5.5
546   79  BH   5.330013    H18B  2023-03-17    H1   8    B    3   9       6.5
547   79  BH   5.330013   H114A  2023-03-17    H1  14    A    5   3       6.2
548   79  BH   5.330013   H117A  2023-03-17    H1  17    A    4   6       3.6
549   79  BH   5.330013   H118B  2023-03-17    H1  18    B    5  18       5.0
550   79  BH   5.330013   H224B  2023-03-17    H2  24    B    5  29       8.5
551   79  BH   5.330013   H233B  2023-03-17    H2  33    B    4   7       3.5
552   79  BH   5.330013   H239B  2023-03-17    H2  39    B    6   2       4.0
553   79  BH   5.330013   H110C  2023-03-17    H1  10    C    6  11       5.6
554   79  BH   5.330013   H115D  2023-03-17    H1  15    D    5  13       5.6
555   79  BH   5.330013   H119C  2023-03-17    H1  19    C    7  21       6.5
556   79  BH   5.330013   H222D  2023-03-17    H2  22    D    2   3       4.4
557   79  BH   5.330013   H223C  2023-03-17    H2  23    C    7  22       4.9
558   79  BH   5.330013   H235C  2023-03-17    H2  35    C    2  13       4.8
559   79  BH   5.330013    J14B  2023-03-17    J1   4    B    2   4       4.5
560   79  BH   5.330013   J119B  2023-03-17    J1  19    B    3  11       5.3
561   79  BH   5.330013   J233A  2023-03-17    J2  33    A    3   2       4.7
562   79  BH   5.330013   J237B  2023-03-17    J2  37    B    1   7       5.3
563   79  BH   5.330013    J19C  2023-03-17    J1   9    C    5  30       8.9
564   79  BH   5.330013   J110C  2023-03-17    J1  10    C    3   1       3.9
565   79  BH   5.330013   J114D  2023-03-17    J1  14    D    4   8       3.8
566   79  BH   5.330013   J119C  2023-03-17    J1  19    C    5  20       7.1
567   79  BH   5.330013   J223D  2023-03-17    J2  23    D    4   9       5.2
568   79  BH   5.330013   J225D  2023-03-17    J2  25    D    2   5       3.0
569   79  BH   5.330013   J230C  2023-03-17    J2  30    C    4  14       8.3
570   79  BH   5.330013   J231C  2023-03-17    J2  31    C    5   6       9.2
571   79  BH   5.330013   J232D  2023-03-17    J2  32    D    3  12       4.5
572   79  BH   5.330013    L15B  2023-03-17    L1   5    B    2   6       4.2
573   79  BH   5.330013   L110A  2023-03-17    L1  10    A    2  14       2.0
574   79  BH   5.330013    L14D  2023-03-17    L1   4    D    4  10       3.8
575   79  BH   5.330013   L113C  2023-03-17    L1  13    C    3  13       3.6
576   79  BH   5.330013   L113D  2023-03-17    L1  13    D    5  22       6.5
577   79  BH   5.330013   L116C  2023-03-17    L1  16    C    6   5       7.9
578   79  BH   5.330013   L228D  2023-03-17    L2  28    D    2   7       5.1
579   79  BH   5.330013   L229D  2023-03-17    L2  29    D    3  14       3.1
580   79  BH   5.330013   L230C  2023-03-17    L2  30    C    4  11       6.6
581   79  BH   5.330013   L234C  2023-03-17    L2  34    C    5   8       4.9
582   79  BH   5.330013   L240D  2023-03-17    L2  40    D    7  13       6.0
583   79  BH   5.330013   L241D  2023-03-17    L2  41    D    4   1       4.6
584   79 DPR   6.412738   F232A  2023-03-17    F2  32    A    6   2       2.5
585   79 DPR   6.412738    F16C  2023-03-17    F1   6    C    6   9       6.6
586   79 DPR   6.412738   F240C  2023-03-17    F2  40    C    7   4      10.2
587   79 DPR   6.412738   H117B  2023-03-17    H1  17    B    5   2      17.0
588   79 DPR   6.412738    H17C  2023-03-17    H1   7    C    2  10      11.4
589   79 DPR   6.412738   H225C  2023-03-17    H2  25    C    2   6      12.4
590   79 DPR   6.412738    J16B  2023-03-17    J1   6    B    6  10       7.4
591   79 DPR   6.412738   J232B  2023-03-17    J2  32    B    6  11       4.2
592   79 DPR   6.412738   L237B  2023-03-17    L2  37    B    6   3       6.6
593   80  CC   5.463058    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
594   80  BH   5.424160    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
595   80 DPR   6.378404    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
596   81  CC   5.936431    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
597   81  BH   5.610675    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
598   81 DPR   6.600084    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
599   82  CC   5.885694    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
600   82  BH   5.678230    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
601   82 DPR   7.017616    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
602   83  CC   5.757422    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
603   83  BH   5.720385    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
604   83 DPR   6.964984    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
605   84  CC   6.028949    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
606   84  BH   5.914924    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
607   84 DPR   7.271102    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
608   85  CC   5.904224    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
609   85  BH   6.021778    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
610   85 DPR   7.429502    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
611   86  CC   6.537206   D112A  2023-03-24    D1  12    A    7   4       5.5
612   86  CC   6.537206   H111B  2023-03-24    H1  11    B    7   2      12.6
613   86  CC   6.537206   H115B  2023-03-24    H1  15    B    3   3      12.5
614   86  CC   6.537206   H116A  2023-03-24    H1  16    A    4  10      10.1
615   86  CC   6.537206   H234B  2023-03-24    H2  34    B    4   7       3.1
616   86  CC   6.537206   H236C  2023-03-24    H2  36    C    2   1       6.0
617   86  CC   6.537206   J112B  2023-03-24    J1  12    B    4   4       9.8
618   86  CC   6.537206   J236C  2023-03-24    J2  36    C    5   5       7.0
619   86  CC   6.537206    L18A  2023-03-24    L1   8    A    4   6       4.6
620   86  CC   6.537206   L232A  2023-03-24    L2  32    A    5   3       6.3
621   86  CC   6.537206   L117C  2023-03-24    L1  17    C    7   3       5.5
622   86  BH   6.080695   D226B  2023-03-24    D2  26    B    5  15       5.4
623   86  BH   6.080695   D230B  2023-03-24    D2  30    B    4   3       3.3
624   86  BH   6.080695   D233B  2023-03-24    D2  33    B    6   8       3.9
625   86  BH   6.080695   D119D  2023-03-24    D1  19    D    5   9       3.5
626   86  BH   6.080695   D229C  2023-03-24    D2  29    C    4  13       4.8
627   86  BH   6.080695   D229D  2023-03-24    D2  29    D    2   9       6.1
628   86  BH   6.080695   D235D  2023-03-24    D2  35    D    3   6       5.7
629   86  BH   6.080695   F115A  2023-03-24    F1  15    A    6   9       4.5
630   86  BH   6.080695   F229B  2023-03-24    F2  29    B    4   5       5.6
631   86  BH   6.080695   F240B  2023-03-24    F2  40    B    5  12       5.7
632   86  BH   6.080695   F115C  2023-03-24    F1  15    C    4   4       4.9
633   86  BH   6.080695   F235D  2023-03-24    F2  35    D    2   1       5.1
634   86  BH   6.080695    H14B  2023-03-24    H1   4    B    2  12       7.3
635   86  BH   6.080695    H18B  2023-03-24    H1   8    B    3   9       5.4
636   86  BH   6.080695   H114A  2023-03-24    H1  14    A    5   3       7.1
637   86  BH   6.080695   H224B  2023-03-24    H2  24    B    5  29       9.8
638   86  BH   6.080695   H233B  2023-03-24    H2  33    B    4   7       4.5
639   86  BH   6.080695   H110C  2023-03-24    H1  10    C    6  11       8.0
640   86  BH   6.080695   H115D  2023-03-24    H1  15    D    5  13       4.9
641   86  BH   6.080695   H119C  2023-03-24    H1  19    C    7  21       7.0
642   86  BH   6.080695   H222D  2023-03-24    H2  22    D    2   3       5.5
643   86  BH   6.080695   H223C  2023-03-24    H2  23    C    7  22       7.4
644   86  BH   6.080695   H235C  2023-03-24    H2  35    C    2  13       3.8
645   86  BH   6.080695    J14B  2023-03-24    J1   4    B    2   4       3.9
646   86  BH   6.080695   J119B  2023-03-24    J1  19    B    3  11       6.4
647   86  BH   6.080695   J233A  2023-03-24    J2  33    A    3   2       5.9
648   86  BH   6.080695   J237B  2023-03-24    J2  37    B    1   7       7.2
649   86  BH   6.080695    J19C  2023-03-24    J1   9    C    5  30       9.7
650   86  BH   6.080695   J110C  2023-03-24    J1  10    C    3   1       4.5
651   86  BH   6.080695   J119C  2023-03-24    J1  19    C    5  20       7.5
652   86  BH   6.080695   J223D  2023-03-24    J2  23    D    4   9       5.5
653   86  BH   6.080695   J225D  2023-03-24    J2  25    D    2   5       3.0
654   86  BH   6.080695   J230C  2023-03-24    J2  30    C    4  14       8.5
655   86  BH   6.080695   J231D  2023-03-24    J2  31    D    7  14       8.3
656   86  BH   6.080695   J232D  2023-03-24    J2  32    D    3  12       4.8
657   86  BH   6.080695    L15B  2023-03-24    L1   5    B    2   6       4.5
658   86  BH   6.080695   L110A  2023-03-24    L1  10    A    2  14       3.1
659   86  BH   6.080695    L14D  2023-03-24    L1   4    D    4  10       4.1
660   86  BH   6.080695   L113C  2023-03-24    L1  13    C    3  13       2.8
661   86  BH   6.080695   L113D  2023-03-24    L1  13    D    5  22       6.5
662   86  BH   6.080695   L228D  2023-03-24    L2  28    D    2   7       5.5
663   86  BH   6.080695   L229D  2023-03-24    L2  29    D    3  14       3.3
664   86  BH   6.080695   L230C  2023-03-24    L2  30    C    4  11       5.1
665   86  BH   6.080695   L240D  2023-03-24    L2  40    D    7  13       5.3
666   86  BH   6.080695   L241D  2023-03-24    L2  41    D    4   1       4.7
667   86 DPR   7.874193    F16C  2023-03-24    F1   6    C    6   9       9.1
668   86 DPR   7.874193   F240C  2023-03-24    F2  40    C    7   4      15.4
669   86 DPR   7.874193   H117B  2023-03-24    H1  17    B    5   2      21.7
670   86 DPR   7.874193    H17C  2023-03-24    H1   7    C    2  10      12.5
671   86 DPR   7.874193   H225C  2023-03-24    H2  25    C    2   6      18.4
672   86 DPR   7.874193    J16B  2023-03-24    J1   6    B    6  10      10.8
673   86 DPR   7.874193   J232B  2023-03-24    J2  32    B    6  11       4.4
674   86 DPR   7.874193   L237B  2023-03-24    L2  37    B    6   3       6.4
675   87  CC   6.837241    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
676   87  BH   6.352095    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
677   87 DPR   8.275922    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
678   88  CC   6.686768    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
679   88  BH   6.237436    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
680   88 DPR   8.207329    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
681   89  CC   7.028307    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
682   89  BH   6.590723    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
683   89 DPR   8.466765    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
684   90  CC   6.955139    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
685   90  BH   6.442395    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
686   90 DPR   8.946546    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
687   91  CC   7.407780    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
688   91  BH   6.252713    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
689   91 DPR   9.472092    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
690   92  CC   7.665776    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
691   92  BH   7.298624    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
692   92 DPR   9.836722    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
693   93  CC   7.930290    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
694   93  BH   7.437636    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
695   93 DPR  10.096792    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
696   94  CC   8.105633    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
697   94  BH   7.557276    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
698   94 DPR  10.153591    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
699   95  CC   8.286193    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
700   95  BH   7.670049    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
701   95 DPR  10.840840    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
702   96  CC   8.789005   H111B  2023-04-03    H1  11    B    7   2      12.1
703   96  CC   8.789005   H115B  2023-04-03    H1  15    B    3   3      14.9
704   96  CC   8.789005   H116A  2023-04-03    H1  16    A    4  10      10.2
705   96  CC   8.789005   J112B  2023-04-03    J1  12    B    4   4       7.5
706   96  CC   8.789005   J236C  2023-04-03    J2  36    C    5   5       8.3
707   96  CC   8.789005   L117C  2023-04-03    L1  17    C    7   3       6.2
708   96  BH   8.164295   D226B  2023-04-03    D2  26    B    5  15       4.9
709   96  BH   8.164295   D230B  2023-04-03    D2  30    B    4   3       3.8
710   96  BH   8.164295   D119D  2023-04-03    D1  19    D    5   9       3.2
711   96  BH   8.164295   D229C  2023-04-03    D2  29    C    4  13       4.9
712   96  BH   8.164295   D229D  2023-04-03    D2  29    D    2   9       4.1
713   96  BH   8.164295   D235D  2023-04-03    D2  35    D    3   6       6.6
714   96  BH   8.164295   F115A  2023-04-03    F1  15    A    6   9       9.9
715   96  BH   8.164295   F229B  2023-04-03    F2  29    B    4   5       5.9
716   96  BH   8.164295   F240B  2023-04-03    F2  40    B    5  12       5.6
717   96  BH   8.164295   F115C  2023-04-03    F1  15    C    4   4       5.9
718   96  BH   8.164295   F235D  2023-04-03    F2  35    D    2   1       5.5
719   96  BH   8.164295    H14B  2023-04-03    H1   4    B    2  12       0.7
720   96  BH   8.164295    H18B  2023-04-03    H1   8    B    3   9       7.1
721   96  BH   8.164295   H114A  2023-04-03    H1  14    A    5   3       8.4
722   96  BH   8.164295   H224B  2023-04-03    H2  24    B    5  29       9.4
723   96  BH   8.164295   H233B  2023-04-03    H2  33    B    4   7       4.8
724   96  BH   8.164295   H110C  2023-04-03    H1  10    C    6  11      10.6
725   96  BH   8.164295   H119C  2023-04-03    H1  19    C    7  21       6.5
726   96  BH   8.164295   H223C  2023-04-03    H2  23    C    7  22       8.4
727   96  BH   8.164295   H235C  2023-04-03    H2  35    C    2  13       3.8
728   96  BH   8.164295    J14B  2023-04-03    J1   4    B    2   4       5.1
729   96  BH   8.164295   J119B  2023-04-03    J1  19    B    3  11       5.9
730   96  BH   8.164295   J233A  2023-04-03    J2  33    A    3   2       5.8
731   96  BH   8.164295   J237B  2023-04-03    J2  37    B    1   7       4.4
732   96  BH   8.164295    J19C  2023-04-03    J1   9    C    5  30      10.2
733   96  BH   8.164295   J110C  2023-04-03    J1  10    C    3   1       5.1
734   96  BH   8.164295   J119C  2023-04-03    J1  19    C    5  20      10.1
735   96  BH   8.164295   J223D  2023-04-03    J2  23    D    4   9       5.7
736   96  BH   8.164295   J225D  2023-04-03    J2  25    D    2   5       6.3
737   96  BH   8.164295   J230C  2023-04-03    J2  30    C    4  14       9.4
738   96  BH   8.164295   J231D  2023-04-03    J2  31    D    7  14      12.1
739   96  BH   8.164295   J232D  2023-04-03    J2  32    D    3  12       5.2
740   96  BH   8.164295    L15B  2023-04-03    L1   5    B    2   6       3.9
741   96  BH   8.164295   L110A  2023-04-03    L1  10    A    2  14       3.6
742   96  BH   8.164295    L14D  2023-04-03    L1   4    D    4  10       7.1
743   96  BH   8.164295   L113C  2023-04-03    L1  13    C    3  13       3.6
744   96  BH   8.164295   L113D  2023-04-03    L1  13    D    5  22       7.3
745   96  BH   8.164295   L228D  2023-04-03    L2  28    D    2   7       7.1
746   96  BH   8.164295   L229D  2023-04-03    L2  29    D    3  14       4.1
747   96  BH   8.164295   L230C  2023-04-03    L2  30    C    4  11       4.3
748   96  BH   8.164295   L240D  2023-04-03    L2  40    D    7  13       6.6
749   96  BH   8.164295   L241D  2023-04-03    L2  41    D    4   1       5.8
750   96 DPR  11.116171    F16C  2023-04-03    F1   6    C    6   9      15.4
751   96 DPR  11.116171   F240C  2023-04-03    F2  40    C    7   4      22.5
752   96 DPR  11.116171   H117B  2023-04-03    H1  17    B    5   2      25.9
753   96 DPR  11.116171    H17C  2023-04-03    H1   7    C    2  10      13.3
754   96 DPR  11.116171   H225C  2023-04-03    H2  25    C    2   6      27.7
755   97  CC   8.781661    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
756   97  BH   8.342990    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
757   97 DPR  11.332256    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
758   98  CC   9.259853    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
759   98  BH   8.660735    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
760   98 DPR  12.027838    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
761   99  CC   9.481174    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
762   99  BH   8.766512    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
763   99 DPR  12.868052    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
764  100  CC   9.856779    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
765  100  BH   9.289722    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
766  100 DPR  12.873822    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
767  101  CC  10.405864    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
768  101  BH   9.374273    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
769  101 DPR  13.650340    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
770  102  CC  10.534922    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
771  102  BH   9.452671    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
772  102 DPR  13.794474    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
773  103  CC  10.609587   H111B  2023-04-10    H1  11    B    7   2      16.2
774  103  CC  10.609587   H115B  2023-04-10    H1  15    B    3   3      23.8
775  103  CC  10.609587   H116A  2023-04-10    H1  16    A    4  10      12.3
776  103  CC  10.609587   J112B  2023-04-10    J1  12    B    4   4      10.1
777  103  CC  10.609587   J236C  2023-04-10    J2  36    C    5   5       5.9
778  103  BH   9.897299   D226B  2023-04-10    D2  26    B    5  15       5.7
779  103  BH   9.897299   D230B  2023-04-10    D2  30    B    4   3       6.8
780  103  BH   9.897299   D119D  2023-04-10    D1  19    D    5   9       3.6
781  103  BH   9.897299   D229C  2023-04-10    D2  29    C    4  13       5.8
782  103  BH   9.897299   D229D  2023-04-10    D2  29    D    2   9       4.3
783  103  BH   9.897299   D235D  2023-04-10    D2  35    D    3   6       8.4
784  103  BH   9.897299   F115A  2023-04-10    F1  15    A    6   9      17.4
785  103  BH   9.897299   F229B  2023-04-10    F2  29    B    4   5       7.2
786  103  BH   9.897299   F240B  2023-04-10    F2  40    B    5  12       6.6
787  103  BH   9.897299   F115C  2023-04-10    F1  15    C    4   4       5.5
788  103  BH   9.897299   F235D  2023-04-10    F2  35    D    2   1       5.7
789  103  BH   9.897299    H14B  2023-04-10    H1   4    B    2  12       6.1
790  103  BH   9.897299    H18B  2023-04-10    H1   8    B    3   9       8.8
791  103  BH   9.897299   H114A  2023-04-10    H1  14    A    5   3      12.3
792  103  BH   9.897299   H224B  2023-04-10    H2  24    B    5  29      13.4
793  103  BH   9.897299   H233B  2023-04-10    H2  33    B    4   7       5.7
794  103  BH   9.897299   H110C  2023-04-10    H1  10    C    6  11      13.2
795  103  BH   9.897299   H119C  2023-04-10    H1  19    C    7  21       5.1
796  103  BH   9.897299   H223C  2023-04-10    H2  23    C    7  22       8.6
797  103  BH   9.897299   H235C  2023-04-10    H2  35    C    2  13       3.9
798  103  BH   9.897299    J14B  2023-04-10    J1   4    B    2   4       6.0
799  103  BH   9.897299   J119B  2023-04-10    J1  19    B    3  11       6.9
800  103  BH   9.897299   J233A  2023-04-10    J2  33    A    3   2       6.4
801  103  BH   9.897299   J237B  2023-04-10    J2  37    B    1   7      12.2
802  103  BH   9.897299    J19C  2023-04-10    J1   9    C    5  30       8.5
803  103  BH   9.897299   J110C  2023-04-10    J1  10    C    3   1       7.0
804  103  BH   9.897299   J119C  2023-04-10    J1  19    C    5  20      11.4
805  103  BH   9.897299   J225D  2023-04-10    J2  25    D    2   5       7.7
806  103  BH   9.897299   J230C  2023-04-10    J2  30    C    4  14       9.0
807  103  BH   9.897299   J231D  2023-04-10    J2  31    D    7  14      10.2
808  103  BH   9.897299   J232D  2023-04-10    J2  32    D    3  12       6.0
809  103  BH   9.897299   L110A  2023-04-10    L1  10    A    2  14       2.1
810  103  BH   9.897299    L14D  2023-04-10    L1   4    D    4  10       4.1
811  103  BH   9.897299   L113C  2023-04-10    L1  13    C    3  13       2.5
812  103  BH   9.897299   L113D  2023-04-10    L1  13    D    5  22       9.0
813  103  BH   9.897299   L228D  2023-04-10    L2  28    D    2   7       6.3
814  103  BH   9.897299   L229D  2023-04-10    L2  29    D    3  14       3.6
815  103  BH   9.897299   L240D  2023-04-10    L2  40    D    7  13       8.4
816  103  BH   9.897299   L241D  2023-04-10    L2  41    D    4   1       4.4
817  103 DPR  14.689341    F16C  2023-04-10    F1   6    C    6   9      16.6
818  103 DPR  14.689341   H117B  2023-04-10    H1  17    B    5   2      38.2
819  103 DPR  14.689341   H225C  2023-04-10    H2  25    C    2   6      37.4
820  103 DPR  14.689341    J16B  2023-04-10    J1   6    B    6  10      23.6
821  104  CC  11.083230    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
822  104  BH  10.124653    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
823  104 DPR  15.281657    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
824  105  CC  11.447500    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
825  105  BH  10.612240    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
826  105 DPR  15.477213    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
827  106  CC  12.361125    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
828  106  BH  11.016104    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
829  106 DPR  16.645268    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
830  107  CC  12.588790    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
831  107  BH  10.914579    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
832  107 DPR  17.170673    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
833  108  CC  12.847926    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
834  108  BH  11.636420    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
835  108 DPR  17.914141    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
836  109  CC  13.749562    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
837  109  BH  11.901993    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
838  109 DPR  18.911071    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
839  110  CC  13.855781   H111B  2023-04-17    H1  11    B    7   2      14.4
840  110  CC  13.855781   H115B  2023-04-17    H1  15    B    3   3      32.6
841  110  CC  13.855781   H116A  2023-04-17    H1  16    A    4  10      15.1
842  110  CC  13.855781   J112B  2023-04-17    J1  12    B    4   4       9.0
843  110  BH  12.318288   D226B  2023-04-17    D2  26    B    5  15       7.9
844  110  BH  12.318288   D230B  2023-04-17    D2  30    B    4   3       8.1
845  110  BH  12.318288   D119D  2023-04-17    D1  19    D    5   9       3.3
846  110  BH  12.318288   D229C  2023-04-17    D2  29    C    4  13       6.8
847  110  BH  12.318288   D229D  2023-04-17    D2  29    D    2   9       7.1
848  110  BH  12.318288   D235D  2023-04-17    D2  35    D    3   6      10.2
849  110  BH  12.318288   F115A  2023-04-17    F1  15    A    6   9      22.9
850  110  BH  12.318288   F229B  2023-04-17    F2  29    B    4   5       8.2
851  110  BH  12.318288   F240B  2023-04-17    F2  40    B    5  12       8.0
852  110  BH  12.318288   F235D  2023-04-17    F2  35    D    2   1       7.3
853  110  BH  12.318288    H14B  2023-04-17    H1   4    B    2  12      11.4
854  110  BH  12.318288    H18B  2023-04-17    H1   8    B    3   9      15.2
855  110  BH  12.318288   H114A  2023-04-17    H1  14    A    5   3      13.1
856  110  BH  12.318288   H224B  2023-04-17    H2  24    B    5  29      19.2
857  110  BH  12.318288   H233B  2023-04-17    H2  33    B    4   7       6.2
858  110  BH  12.318288   H110C  2023-04-17    H1  10    C    6  11      18.3
859  110  BH  12.318288   H223C  2023-04-17    H2  23    C    7  22      11.8
860  110  BH  12.318288    J14B  2023-04-17    J1   4    B    2   4      12.9
861  110  BH  12.318288   J119B  2023-04-17    J1  19    B    3  11      12.0
862  110  BH  12.318288   J233A  2023-04-17    J2  33    A    3   2       9.0
863  110  BH  12.318288   J237B  2023-04-17    J2  37    B    1   7      18.4
864  110  BH  12.318288    J19C  2023-04-17    J1   9    C    5  30       9.8
865  110  BH  12.318288   J110C  2023-04-17    J1  10    C    3   1      13.4
866  110  BH  12.318288   J119C  2023-04-17    J1  19    C    5  20      18.4
867  110  BH  12.318288   J223D  2023-04-17    J2  23    D    4   9       6.1
868  110  BH  12.318288   J225D  2023-04-17    J2  25    D    2   5       8.8
869  110  BH  12.318288   J230C  2023-04-17    J2  30    C    4  14      13.1
870  110  BH  12.318288   J231D  2023-04-17    J2  31    D    7  14      11.2
871  110  BH  12.318288   J232D  2023-04-17    J2  32    D    3  12      10.8
872  110  BH  12.318288   L110A  2023-04-17    L1  10    A    2  14      11.1
873  110  BH  12.318288    L14D  2023-04-17    L1   4    D    4  10      14.4
874  110  BH  12.318288   L113C  2023-04-17    L1  13    C    3  13      16.0
875  110  BH  12.318288   L113D  2023-04-17    L1  13    D    5  22      12.4
876  110  BH  12.318288   L228D  2023-04-17    L2  28    D    2   7       7.1
877  110  BH  12.318288   L229D  2023-04-17    L2  29    D    3  14       5.6
878  110  BH  12.318288   L240D  2023-04-17    L2  40    D    7  13      11.8
879  110  BH  12.318288   L241D  2023-04-17    L2  41    D    4   1       6.9
880  110 DPR  19.440602    F16C  2023-04-17    F1   6    C    6   9      16.1
881  110 DPR  19.440602   H117B  2023-04-17    H1  17    B    5   2      42.7
882  110 DPR  19.440602   H225C  2023-04-17    H2  25    C    2   6      36.2
883  110 DPR  19.440602    J16B  2023-04-17    J1   6    B    6  10      24.5
884  111  CC  14.596189    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
885  111  BH  13.114533    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
886  111 DPR  20.102806    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
887  112  CC  14.984914    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
888  112  BH  13.349971    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
889  112 DPR  21.234189    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
890  113  CC  15.368785    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
891  113  BH  13.668082    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
892  113 DPR  21.784671    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
893  114  CC  16.574341    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
894  114  BH  14.612022    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
895  114 DPR  22.674816    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
896  115  CC  17.282482    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
897  115  BH  14.791267    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
898  115 DPR  23.194299    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
899  116  CC  17.065289    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
900  116  BH  15.272846    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
901  116 DPR  24.656664    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
902  117  CC  17.878895   H111B  2023-04-24    H1  11    B    7   2      27.1
903  117  CC  17.878895   H115B  2023-04-24    H1  15    B    3   3      37.0
904  117  CC  17.878895   H116A  2023-04-24    H1  16    A    4  10      15.9
905  117  BH  15.991544   D226B  2023-04-24    D2  26    B    5  15      15.9
906  117  BH  15.991544   D230B  2023-04-24    D2  30    B    4   3      12.6
907  117  BH  15.991544   D119D  2023-04-24    D1  19    D    5   9       2.1
908  117  BH  15.991544   D229C  2023-04-24    D2  29    C    4  13       8.1
909  117  BH  15.991544   D229D  2023-04-24    D2  29    D    2   9      10.9
910  117  BH  15.991544   D235D  2023-04-24    D2  35    D    3   6      16.3
911  117  BH  15.991544   F115A  2023-04-24    F1  15    A    6   9      33.1
912  117  BH  15.991544   F229B  2023-04-24    F2  29    B    4   5       9.1
913  117  BH  15.991544   F240B  2023-04-24    F2  40    B    5  12      10.1
914  117  BH  15.991544   F235D  2023-04-24    F2  35    D    2   1       1.8
915  117  BH  15.991544    H14B  2023-04-24    H1   4    B    2  12      14.9
916  117  BH  15.991544    H18B  2023-04-24    H1   8    B    3   9      21.4
917  117  BH  15.991544   H114A  2023-04-24    H1  14    A    5   3      17.5
918  117  BH  15.991544   H224B  2023-04-24    H2  24    B    5  29      28.3
919  117  BH  15.991544   H233B  2023-04-24    H2  33    B    4   7      13.7
920  117  BH  15.991544   H110C  2023-04-24    H1  10    C    6  11      24.4
921  117  BH  15.991544   H223C  2023-04-24    H2  23    C    7  22      12.2
922  117  BH  15.991544    J14B  2023-04-24    J1   4    B    2   4       5.2
923  117  BH  15.991544   J119B  2023-04-24    J1  19    B    3  11      11.9
924  117  BH  15.991544   J233A  2023-04-24    J2  33    A    3   2       7.9
925  117  BH  15.991544   J237B  2023-04-24    J2  37    B    1   7      25.3
926  117  BH  15.991544   J110C  2023-04-24    J1  10    C    3   1      14.8
927  117  BH  15.991544   J119C  2023-04-24    J1  19    C    5  20      28.3
928  117  BH  15.991544   J223D  2023-04-24    J2  23    D    4   9       8.4
929  117  BH  15.991544   J225D  2023-04-24    J2  25    D    2   5      14.1
930  117  BH  15.991544   J230C  2023-04-24    J2  30    C    4  14      15.3
931  117  BH  15.991544   J231D  2023-04-24    J2  31    D    7  14      10.6
932  117  BH  15.991544   J232D  2023-04-24    J2  32    D    3  12      19.2
933  117  BH  15.991544   L110A  2023-04-24    L1  10    A    2  14      14.4
934  117  BH  15.991544    L14D  2023-04-24    L1   4    D    4  10      20.6
935  117  BH  15.991544   L113C  2023-04-24    L1  13    C    3  13      27.7
936  117  BH  15.991544   L113D  2023-04-24    L1  13    D    5  22      17.5
937  117  BH  15.991544   L228D  2023-04-24    L2  28    D    2   7       6.9
938  117  BH  15.991544   L229D  2023-04-24    L2  29    D    3  14      11.8
939  117  BH  15.991544   L240D  2023-04-24    L2  40    D    7  13      14.1
940  117  BH  15.991544   L241D  2023-04-24    L2  41    D    4   1      12.3
941  117 DPR  25.456961    F16C  2023-04-24    F1   6    C    6   9      16.0
942  117 DPR  25.456961   H117B  2023-04-24    H1  17    B    5   2      42.3
943  117 DPR  25.456961   H225C  2023-04-24    H2  25    C    2   6      36.1
944  117 DPR  25.456961    J16B  2023-04-24    J1   6    B    6  10      26.7
945  118  CC  18.598172    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
946  118  BH  16.563027    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
947  118 DPR  26.663001    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
948  119  CC  19.364915    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
949  119  BH  17.142898    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
950  119 DPR  27.538988    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
951  120  CC  19.996646    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
952  120  BH  17.764244    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
953  120 DPR  28.327148    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
954  121  CC  20.951338    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
955  121  BH  18.036483    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
956  121 DPR  29.721795    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
957  122  CC  21.328928    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
958  122  BH  18.949765    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
959  122 DPR  31.203732    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
960  123  CC  22.019455    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
961  123  BH  19.625868    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
962  123 DPR  31.904276    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
963  124  CC  22.980478   H111B  2023-05-01    H1  11    B    7   2      15.5
964  124  CC  22.980478   H115B  2023-05-01    H1  15    B    3   3      34.7
965  124  CC  22.980478   H116A  2023-05-01    H1  16    A    4  10      13.5
966  124  BH  20.192369   D226B  2023-05-01    D2  26    B    5  15      20.5
967  124  BH  20.192369   D230B  2023-05-01    D2  30    B    4   3      17.2
968  124  BH  20.192369   D119D  2023-05-01    D1  19    D    5   9       3.1
969  124  BH  20.192369   D229C  2023-05-01    D2  29    C    4  13      10.8
970  124  BH  20.192369   D229D  2023-05-01    D2  29    D    2   9      13.1
971  124  BH  20.192369   D235D  2023-05-01    D2  35    D    3   6       3.9
972  124  BH  20.192369   F115A  2023-05-01    F1  15    A    6   9      42.5
973  124  BH  20.192369   F229B  2023-05-01    F2  29    B    4   5      12.5
974  124  BH  20.192369   F240B  2023-05-01    F2  40    B    5  12      16.5
975  124  BH  20.192369   F235D  2023-05-01    F2  35    D    2   1      10.2
976  124  BH  20.192369    H14B  2023-05-01    H1   4    B    2  12      18.6
977  124  BH  20.192369    H18B  2023-05-01    H1   8    B    3   9      30.5
978  124  BH  20.192369   H114A  2023-05-01    H1  14    A    5   3      19.5
979  124  BH  20.192369   H224B  2023-05-01    H2  24    B    5  29      38.4
980  124  BH  20.192369   H233B  2023-05-01    H2  33    B    4   7      19.1
981  124  BH  20.192369   H110C  2023-05-01    H1  10    C    6  11      29.9
982  124  BH  20.192369   H223C  2023-05-01    H2  23    C    7  22      12.1
983  124  BH  20.192369    J14B  2023-05-01    J1   4    B    2   4      11.5
984  124  BH  20.192369   J119B  2023-05-01    J1  19    B    3  11      18.4
985  124  BH  20.192369   J233A  2023-05-01    J2  33    A    3   2       9.5
986  124  BH  20.192369   J237B  2023-05-01    J2  37    B    1   7      33.0
987  124  BH  20.192369   J110C  2023-05-01    J1  10    C    3   1      31.5
988  124  BH  20.192369   J119C  2023-05-01    J1  19    C    5  20      44.3
989  124  BH  20.192369   J223D  2023-05-01    J2  23    D    4   9      10.2
990  124  BH  20.192369   J225D  2023-05-01    J2  25    D    2   5      17.0
991  124  BH  20.192369   J230C  2023-05-01    J2  30    C    4  14      23.5
992  124  BH  20.192369   J231D  2023-05-01    J2  31    D    7  14      10.5
993  124  BH  20.192369   J232D  2023-05-01    J2  32    D    3  12      23.5
994  124  BH  20.192369   L110A  2023-05-01    L1  10    A    2  14      21.0
995  124  BH  20.192369    L14D  2023-05-01    L1   4    D    4  10      24.7
996  124  BH  20.192369   L113C  2023-05-01    L1  13    C    3  13      50.6
997  124  BH  20.192369   L113D  2023-05-01    L1  13    D    5  22      21.1
998  124  BH  20.192369   L228D  2023-05-01    L2  28    D    2   7       7.5
999  124  BH  20.192369   L229D  2023-05-01    L2  29    D    3  14      16.5
1000 124  BH  20.192369   L240D  2023-05-01    L2  40    D    7  13      15.0
1001 124  BH  20.192369   L241D  2023-05-01    L2  41    D    4   1      16.5
1002 124 DPR  33.059192    F16C  2023-05-01    F1   6    C    6   9      17.0
1003 124 DPR  33.059192   H117B  2023-05-01    H1  17    B    5   2      42.4
1004 124 DPR  33.059192    J16B  2023-05-01    J1   6    B    6  10      23.5
1005 125  CC  23.829158    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1006 125  BH  20.901241    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1007 125 DPR  34.563889    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1008 126  CC  24.562509    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1009 126  BH  21.545037    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1010 126 DPR  35.559782    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1011 127  CC  25.882807    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1012 127  BH  22.143322    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1013 127 DPR  36.798499    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1014 128  CC  26.097965    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1015 128  BH  23.121734    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1016 128 DPR  37.936337    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1017 129  CC  26.943909    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1018 129  BH  23.543188    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1019 129 DPR  39.534618    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1020 130  CC  28.060317    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1021 130  BH  24.673691    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1022 130 DPR  41.134888    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1023 131  CC  29.187596   H116A  2023-05-08    H1  16    A    4  10      15.1
1024 131  BH  25.238720   D226B  2023-05-08    D2  26    B    5  15      23.1
1025 131  BH  25.238720   D230B  2023-05-08    D2  30    B    4   3      22.1
1026 131  BH  25.238720   D119D  2023-05-08    D1  19    D    5   9       3.1
1027 131  BH  25.238720   D229C  2023-05-08    D2  29    C    4  13      13.3
1028 131  BH  25.238720   D229D  2023-05-08    D2  29    D    2   9      15.8
1029 131  BH  25.238720   D235D  2023-05-08    D2  35    D    3   6      35.2
1030 131  BH  25.238720   F115A  2023-05-08    F1  15    A    6   9      51.3
1031 131  BH  25.238720   F229B  2023-05-08    F2  29    B    4   5      14.6
1032 131  BH  25.238720   F240B  2023-05-08    F2  40    B    5  12      18.7
1033 131  BH  25.238720   F235D  2023-05-08    F2  35    D    2   1      12.2
1034 131  BH  25.238720    H14B  2023-05-08    H1   4    B    2  12      20.8
1035 131  BH  25.238720    H18B  2023-05-08    H1   8    B    3   9      36.1
1036 131  BH  25.238720   H114A  2023-05-08    H1  14    A    5   3      21.6
1037 131  BH  25.238720   H224B  2023-05-08    H2  24    B    5  29      48.3
1038 131  BH  25.238720   H233B  2023-05-08    H2  33    B    4   7      20.8
1039 131  BH  25.238720   H110C  2023-05-08    H1  10    C    6  11      34.0
1040 131  BH  25.238720   H223C  2023-05-08    H2  23    C    7  22      12.1
1041 131  BH  25.238720    J14B  2023-05-08    J1   4    B    2   4       4.4
1042 131  BH  25.238720   J119B  2023-05-08    J1  19    B    3  11      18.6
1043 131  BH  25.238720   J233A  2023-05-08    J2  33    A    3   2       5.7
1044 131  BH  25.238720   J237B  2023-05-08    J2  37    B    1   7      39.2
1045 131  BH  25.238720   J110C  2023-05-08    J1  10    C    3   1      33.1
1046 131  BH  25.238720   J119C  2023-05-08    J1  19    C    5  20      60.5
1047 131  BH  25.238720   J223D  2023-05-08    J2  23    D    4   9      13.2
1048 131  BH  25.238720   J225D  2023-05-08    J2  25    D    2   5      20.4
1049 131  BH  25.238720   J230C  2023-05-08    J2  30    C    4  14      24.6
1050 131  BH  25.238720   J231D  2023-05-08    J2  31    D    7  14       8.9
1051 131  BH  25.238720   J232D  2023-05-08    J2  32    D    3  12      24.0
1052 131  BH  25.238720   L110A  2023-05-08    L1  10    A    2  14      23.3
1053 131  BH  25.238720    L14D  2023-05-08    L1   4    D    4  10      27.3
1054 131  BH  25.238720   L113C  2023-05-08    L1  13    C    3  13      65.9
1055 131  BH  25.238720   L113D  2023-05-08    L1  13    D    5  22      21.2
1056 131  BH  25.238720   L228D  2023-05-08    L2  28    D    2   7       5.0
1057 131  BH  25.238720   L229D  2023-05-08    L2  29    D    3  14      17.8
1058 131  BH  25.238720   L240D  2023-05-08    L2  40    D    7  13      14.1
1059 131  BH  25.238720   L241D  2023-05-08    L2  41    D    4   1      11.5
1060 131 DPR  42.506376   H117B  2023-05-08    H1  17    B    5   2      42.2
1061 132  CC  30.298348    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1062 132  BH  26.397878    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1063 132 DPR  43.996632    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1064 133  CC  30.952394    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1065 133  BH  26.882530    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1066 133 DPR  45.612730    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1067 134  CC  32.107878    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1068 134  BH  27.963245    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1069 134 DPR  46.709056    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1070 135  CC  32.877267    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1071 135  BH  28.593949    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1072 135 DPR  47.995706    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1073 136  CC  33.739528    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1074 136  BH  29.542303    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1075 136 DPR  49.916577    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1076 137  CC  34.923623    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1077 137  BH  30.122098    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1078 137 DPR  51.429858    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1079 138  CC  35.660916    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1080 138  BH  31.196302    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1081 138 DPR  52.695856    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1082 139  CC  36.823764    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1083 139  BH  32.127654    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1084 139 DPR  54.657727    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1085 140  CC  37.850210    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1086 140  BH  33.075873    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1087 140 DPR  56.026651    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1088 141  CC  38.871460    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1089 141  BH  34.197833    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1090 141 DPR  57.420682    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1091 142  CC  39.953462    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1092 142  BH  34.706582    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1093 142 DPR  59.065748    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1094 143  CC  41.037084    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1095 143  BH  35.568389    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1096 143 DPR  60.919146    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1097 144  CC  42.137516    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1098 144  BH  36.389452    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1099 144 DPR  62.641743    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1100 145  CC  43.366689    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1101 145  BH  37.669593    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1102 145 DPR  64.242879    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1103 146  CC  44.541872    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1104 146  BH  38.567430    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1105 146 DPR  65.630037    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1106 147  CC  45.444569    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1107 147  BH  39.422887    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1108 147 DPR  67.617845    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1109 148  CC  46.666651    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1110 148  BH  40.206033    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1111 148 DPR  68.903683    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1112 149  CC  47.288076    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1113 149  BH  41.310744    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1114 149 DPR  70.554517    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1115 150  CC  48.363093    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1116 150  BH  42.133949    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1117 150 DPR  72.445219    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1118 151  CC  49.780894    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1119 151  BH  43.130438    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1120 151 DPR  73.969456    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1121 152  CC  50.666027    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1122 152  BH  43.666691    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1123 152 DPR  75.740001    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1124 153  CC  51.951086    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1125 153  BH  44.869288    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1126 153 DPR  77.005471    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1127 154  CC  52.510411    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1128 154  BH  45.680591    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1129 154 DPR  78.770275    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1130 155  CC  54.278327    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1131 155  BH  46.380003    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1132 155 DPR  80.366905    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1133 156  CC  54.625418    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1134 156  BH  47.619006    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1135 156 DPR  81.784463    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1136 157  CC  55.572277    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1137 157  BH  48.121739    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1138 157 DPR  83.408759    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1139 158  CC  56.972035    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1140 158  BH  49.071728    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1141 158 DPR  84.916182    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1142 159  CC  57.733157    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1143 159  BH  50.085899   D226B  2023-06-05    D2  26    B    5  15      77.4
1144 159  BH  50.085899   D230B  2023-06-05    D2  30    B    4   3      52.1
1145 159  BH  50.085899   D229C  2023-06-05    D2  29    C    4  13      23.7
1146 159  BH  50.085899   D229D  2023-06-05    D2  29    D    2   9      34.5
1147 159  BH  50.085899   D235D  2023-06-05    D2  35    D    3   6      80.5
1148 159  BH  50.085899   F115A  2023-06-05    F1  15    A    6   9      88.9
1149 159  BH  50.085899   F229B  2023-06-05    F2  29    B    4   5      26.3
1150 159  BH  50.085899   F240B  2023-06-05    F2  40    B    5  12      47.0
1151 159  BH  50.085899   F235D  2023-06-05    F2  35    D    2   1      37.5
1152 159  BH  50.085899    H14B  2023-06-05    H1   4    B    2  12      28.2
1153 159  BH  50.085899    H18B  2023-06-05    H1   8    B    3   9      76.1
1154 159  BH  50.085899   H114A  2023-06-05    H1  14    A    5   3      19.9
1155 159  BH  50.085899   H224B  2023-06-05    H2  24    B    5  29      70.8
1156 159  BH  50.085899   H233B  2023-06-05    H2  33    B    4   7      21.1
1157 159  BH  50.085899   H110C  2023-06-05    H1  10    C    6  11      44.8
1158 159  BH  50.085899   J119B  2023-06-05    J1  19    B    3  11      31.9
1159 159  BH  50.085899   J237B  2023-06-05    J2  37    B    1   7      39.9
1160 159  BH  50.085899   J110C  2023-06-05    J1  10    C    3   1      70.0
1161 159  BH  50.085899   J119C  2023-06-05    J1  19    C    5  20     102.4
1162 159  BH  50.085899   J223D  2023-06-05    J2  23    D    4   9      51.4
1163 159  BH  50.085899   J225D  2023-06-05    J2  25    D    2   5      71.0
1164 159  BH  50.085899   J230C  2023-06-05    J2  30    C    4  14      32.3
1165 159  BH  50.085899   J232D  2023-06-05    J2  32    D    3  12      27.8
1166 159  BH  50.085899   L110A  2023-06-05    L1  10    A    2  14      61.7
1167 159  BH  50.085899    L14D  2023-06-05    L1   4    D    4  10      77.6
1168 159  BH  50.085899   L113C  2023-06-05    L1  13    C    3  13     120.6
1169 159  BH  50.085899   L113D  2023-06-05    L1  13    D    5  22      58.9
1170 159  BH  50.085899   L229D  2023-06-05    L2  29    D    3  14      31.4
1171 159  BH  50.085899   L240D  2023-06-05    L2  40    D    7  13      19.7
1172 159  BH  50.085899   L241D  2023-06-05    L2  41    D    4   1      49.4
1173 159 DPR  85.902175    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1174 160  CC  58.310997    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1175 160  BH  50.620510    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1176 160 DPR  87.692400    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1177 161  CC  59.530560    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1178 161  BH  51.393821    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1179 161 DPR  89.118147    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1180 162  CC  60.792405    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1181 162  BH  52.307221    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1182 162 DPR  90.562057    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1183 163  CC  61.234103    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1184 163  BH  53.149802    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1185 163 DPR  92.080028    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1186 164  CC  62.209131    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1187 164  BH  53.733270    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1188 164 DPR  93.121142    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1189 165  CC  62.832858    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1190 165  BH  54.434346    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1191 165 DPR  94.418714    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1192 166  CC  63.637564    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1193 166  BH  55.233822    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1194 166 DPR  95.440624    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1195 167  CC  64.466386    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1196 167  BH  55.805459    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1197 167 DPR  96.809681    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1198 168  CC  65.057894    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1199 168  BH  56.316710    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1200 168 DPR  97.977634    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1201 169  CC  65.782181    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1202 169  BH  56.924255    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1203 169 DPR  98.605446    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1204 170  CC  66.296381    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1205 170  BH  57.311219    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1206 170 DPR 100.135224    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1207 171  CC  67.136209    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1208 171  BH  58.141660    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1209 171 DPR 101.072480    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1210 172  CC  67.922378    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1211 172  BH  58.303948    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1212 172 DPR 102.289417    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1213 173  CC  68.715501    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1214 173  BH  59.313891    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1215 173 DPR 102.600015    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1216 174  CC  69.016414    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1217 174  BH  59.459249    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1218 174 DPR 104.090764    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1219 175  CC  69.336516    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1220 175  BH  60.121786    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1221 175 DPR 104.873503    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1222 176  CC  69.800120    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1223 176  BH  60.617953    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1224 176 DPR 105.364535    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1225 177  CC  70.478293    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1226 177  BH  60.924538    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1227 177 DPR 106.117513    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1228 178  CC  70.938890    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1229 178  BH  60.904400    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1230 178 DPR 106.831206    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1231 179  CC  71.779977    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1232 179  BH  61.725652    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1233 179 DPR 107.007783    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1234 180  CC  71.879739    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1235 180  BH  61.509736    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1236 180 DPR 108.176172    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1237 181  CC  71.923660    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1238 181  BH  62.113706    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1239 181 DPR 108.647118    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1240 182  CC  72.356067    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1241 182  BH  62.707041    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1242 182 DPR 109.009752    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1243 183  CC  72.627308    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1244 183  BH  62.739367    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1245 183 DPR 109.543316    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1246 184  CC  72.992116    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1247 184  BH  63.047961    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1248 184 DPR 109.785786    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1249 185  CC  73.245853    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1250 185  BH  63.482200    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1251 185 DPR 110.206729    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1252 186  CC  73.570108    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1253 186  BH  63.409155    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1254 186 DPR 111.136038    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1255 187  CC  73.750035    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1256 187  BH  63.517974   D226B  2023-07-03    D2  26    B    5  15      67.4
1257 187  BH  63.517974   D230B  2023-07-03    D2  30    B    4   3      76.3
1258 187  BH  63.517974   D229C  2023-07-03    D2  29    C    4  13      26.9
1259 187  BH  63.517974   D229D  2023-07-03    D2  29    D    2   9      42.7
1260 187  BH  63.517974   D235D  2023-07-03    D2  35    D    3   6      83.2
1261 187  BH  63.517974   F115A  2023-07-03    F1  15    A    6   9     110.4
1262 187  BH  63.517974   F229B  2023-07-03    F2  29    B    4   5      34.3
1263 187  BH  63.517974   F240B  2023-07-03    F2  40    B    5  12      70.9
1264 187  BH  63.517974   F235D  2023-07-03    F2  35    D    2   1      67.8
1265 187  BH  63.517974    H18B  2023-07-03    H1   8    B    3   9      81.6
1266 187  BH  63.517974   H114A  2023-07-03    H1  14    A    5   3      18.7
1267 187  BH  63.517974   H224B  2023-07-03    H2  24    B    5  29      69.0
1268 187  BH  63.517974   H233B  2023-07-03    H2  33    B    4   7      18.0
1269 187  BH  63.517974   H110C  2023-07-03    H1  10    C    6  11      45.1
1270 187  BH  63.517974   J119B  2023-07-03    J1  19    B    3  11      32.1
1271 187  BH  63.517974   J110C  2023-07-03    J1  10    C    3   1      67.0
1272 187  BH  63.517974   J119C  2023-07-03    J1  19    C    5  20     101.0
1273 187  BH  63.517974   J223D  2023-07-03    J2  23    D    4   9      49.0
1274 187  BH  63.517974   J225D  2023-07-03    J2  25    D    2   5      82.0
1275 187  BH  63.517974   J230C  2023-07-03    J2  30    C    4  14      35.0
1276 187  BH  63.517974   J232D  2023-07-03    J2  32    D    3  12      29.0
1277 187  BH  63.517974   L110A  2023-07-03    L1  10    A    2  14      89.5
1278 187  BH  63.517974    L14D  2023-07-03    L1   4    D    4  10      92.5
1279 187  BH  63.517974   L113C  2023-07-03    L1  13    C    3  13      82.0
1280 187  BH  63.517974   L113D  2023-07-03    L1  13    D    5  22     113.0
1281 187  BH  63.517974   L229D  2023-07-03    L2  29    D    3  14      28.0
1282 187  BH  63.517974   L240D  2023-07-03    L2  40    D    7  13      17.0
1283 187  BH  63.517974   L241D  2023-07-03    L2  41    D    4   1      67.5
1284 187 DPR 111.376556    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1285 188  CC  74.037239    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1286 188  BH  63.996573    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1287 188 DPR 111.119063    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1288 189  CC  74.372764    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1289 189  BH  63.987599    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1290 189 DPR 111.677563    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1291 190  CC  74.057509    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1292 190  BH  64.356749    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1293 190 DPR 112.036224    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1294 191  CC  74.627057    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1295 191  BH  64.100148    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1296 191 DPR 112.562118    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1297 192  CC  74.792647    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1298 192  BH  64.309322    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1299 192 DPR 112.373866    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1300 193  CC  75.150504    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1301 193  BH  64.598085    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1302 193 DPR 112.545399    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1303 194  CC  74.957887    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1304 194  BH  64.645482   D226B  2023-07-10    D2  26    B    5  15      70.0
1305 194  BH  64.645482   D230B  2023-07-10    D2  30    B    4   3      73.0
1306 194  BH  64.645482   D229C  2023-07-10    D2  29    C    4  13      26.5
1307 194  BH  64.645482   D229D  2023-07-10    D2  29    D    2   9      43.0
1308 194  BH  64.645482   D235D  2023-07-10    D2  35    D    3   6      85.8
1309 194  BH  64.645482   F115A  2023-07-10    F1  15    A    6   9     113.2
1310 194  BH  64.645482   F229B  2023-07-10    F2  29    B    4   5      33.2
1311 194  BH  64.645482   F240B  2023-07-10    F2  40    B    5  12      71.6
1312 194  BH  64.645482   F235D  2023-07-10    F2  35    D    2   1      68.1
1313 194  BH  64.645482    H18B  2023-07-10    H1   8    B    3   9      82.2
1314 194  BH  64.645482   H224B  2023-07-10    H2  24    B    5  29      68.1
1315 194  BH  64.645482   H233B  2023-07-10    H2  33    B    4   7      23.5
1316 194  BH  64.645482   J223D  2023-07-10    J2  23    D    4   9      49.6
1317 194  BH  64.645482   J225D  2023-07-10    J2  25    D    2   5      77.0
1318 194  BH  64.645482   J230C  2023-07-10    J2  30    C    4  14      35.1
1319 194  BH  64.645482   J232D  2023-07-10    J2  32    D    3  12      30.5
1320 194  BH  64.645482   L110A  2023-07-10    L1  10    A    2  14      89.4
1321 194  BH  64.645482    L14D  2023-07-10    L1   4    D    4  10      95.6
1322 194  BH  64.645482   L113C  2023-07-10    L1  13    C    3  13     130.0
1323 194  BH  64.645482   L113D  2023-07-10    L1  13    D    5  22      96.1
1324 194  BH  64.645482   L240D  2023-07-10    L2  40    D    7  13      20.1
1325 194  BH  64.645482   L241D  2023-07-10    L2  41    D    4   1      68.4
1326 194 DPR 112.932658    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1327 195  CC  74.735035    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1328 195  BH  64.486604    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1329 195 DPR 113.043086    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1330 196  CC  75.123928    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1331 196  BH  65.006054    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1332 196 DPR 113.275804    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1333 197  CC  74.950744    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1334 197  BH  65.007157    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1335 197 DPR 113.409906    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1336 198  CC  75.390754    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1337 198  BH  65.004296    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1338 198 DPR 113.589558    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1339 199  CC  75.384622    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1340 199  BH  64.960507    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1341 199 DPR 113.362143    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1342 200  CC  75.463249    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1343 200  BH  64.914919    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1344 200 DPR 113.933741    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1345 201  CC  75.795368    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1346 201  BH  65.156515   D226B  2023-07-17    D2  26    B    5  15      68.2
1347 201  BH  65.156515   D230B  2023-07-17    D2  30    B    4   3      77.8
1348 201  BH  65.156515   D229C  2023-07-17    D2  29    C    4  13      25.7
1349 201  BH  65.156515   D229D  2023-07-17    D2  29    D    2   9      44.2
1350 201  BH  65.156515   D235D  2023-07-17    D2  35    D    3   6      83.6
1351 201  BH  65.156515   F115A  2023-07-17    F1  15    A    6   9     113.2
1352 201  BH  65.156515   F229B  2023-07-17    F2  29    B    4   5      34.3
1353 201  BH  65.156515   F240B  2023-07-17    F2  40    B    5  12      76.3
1354 201  BH  65.156515   F235D  2023-07-17    F2  35    D    2   1      68.8
1355 201  BH  65.156515    H18B  2023-07-17    H1   8    B    3   9      83.7
1356 201  BH  65.156515   H224B  2023-07-17    H2  24    B    5  29      70.3
1357 201  BH  65.156515   H233B  2023-07-17    H2  33    B    4   7      22.9
1358 201  BH  65.156515   J223D  2023-07-17    J2  23    D    4   9      50.8
1359 201  BH  65.156515   J225D  2023-07-17    J2  25    D    2   5      83.5
1360 201  BH  65.156515   J230C  2023-07-17    J2  30    C    4  14      41.6
1361 201  BH  65.156515   J232D  2023-07-17    J2  32    D    3  12      30.7
1362 201  BH  65.156515   L110A  2023-07-17    L1  10    A    2  14      90.0
1363 201  BH  65.156515    L14D  2023-07-17    L1   4    D    4  10      73.2
1364 201  BH  65.156515   L113C  2023-07-17    L1  13    C    3  13     107.6
1365 201  BH  65.156515   L241D  2023-07-17    L2  41    D    4   1      37.4
1366 201 DPR 113.922013    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1367 202  CC  75.306502    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1368 202  BH  65.510875    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1369 202 DPR 113.867763    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1370 203  CC  75.628981    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1371 203  BH  64.918779    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1372 203 DPR 114.115781    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1373 204  CC  75.640226    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1374 204  BH  65.274676    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1375 204 DPR 114.114973    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1376 205  CC  75.511065    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1377 205  BH  65.436566    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1378 205 DPR 113.648827    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1379 206  CC  75.701962    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1380 206  BH  65.552906    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1381 206 DPR 114.070903    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1382 207  CC  76.056240    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1383 207  BH  65.800849    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1384 207 DPR 114.395974    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1385 208  CC  76.029601    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
1386 208  BH  65.129865   D226B  2023-07-24    D2  26    B    5  15      68.8
1387 208  BH  65.129865   D230B  2023-07-24    D2  30    B    4   3      79.3
1388 208  BH  65.129865   D229C  2023-07-24    D2  29    C    4  13      26.4
1389 208  BH  65.129865   D229D  2023-07-24    D2  29    D    2   9      44.8
1390 208  BH  65.129865   D235D  2023-07-24    D2  35    D    3   6      86.0
1391 208  BH  65.129865   F115A  2023-07-24    F1  15    A    6   9     109.9
1392 208  BH  65.129865   F229B  2023-07-24    F2  29    B    4   5      36.2
1393 208  BH  65.129865   F240B  2023-07-24    F2  40    B    5  12      75.1
1394 208  BH  65.129865   F235D  2023-07-24    F2  35    D    2   1      69.3
1395 208  BH  65.129865    H18B  2023-07-24    H1   8    B    3   9      81.7
1396 208  BH  65.129865   H224B  2023-07-24    H2  24    B    5  29      68.5
1397 208  BH  65.129865   H233B  2023-07-24    H2  33    B    4   7      22.8
1398 208  BH  65.129865   J225D  2023-07-24    J2  25    D    2   5      81.5
1399 208  BH  65.129865   J230C  2023-07-24    J2  30    C    4  14      32.8
1400 208  BH  65.129865   J232D  2023-07-24    J2  32    D    3  12      30.2
1401 208  BH  65.129865   L110A  2023-07-24    L1  10    A    2  14      91.2
1402 208  BH  65.129865    L14D  2023-07-24    L1   4    D    4  10      93.3
1403 208  BH  65.129865   L113C  2023-07-24    L1  13    C    3  13      89.8
1404 208  BH  65.129865   L113D  2023-07-24    L1  13    D    5  22     112.1
1405 208  BH  65.129865   L241D  2023-07-24    L2  41    D    4   1      69.6
1406 208 DPR 113.929124    <NA>        <NA>  <NA>  NA <NA> <NA>  NA        NA
     herbivory     pred2
1         <NA>  4.205860
2         <NA>  3.796339
3         <NA>  3.926453
4            Y  3.998007
5         <NA>  4.301834
6            Y  3.983936
7            Y  4.078243
8         <NA>  4.156081
9         <NA>  3.968683
10        <NA>  4.166891
11           Y  4.369282
12        <NA>  4.022264
13           Y  4.222716
14        <NA>  3.955638
15        <NA>  3.939511
16        <NA>  4.026712
17           Y  4.104549
18        <NA>  4.068152
19           Y  3.871207
20           Y  4.228728
21        <NA>  4.066312
22        <NA>  3.920191
23        <NA>  3.575278
24        <NA>  4.444030
25        <NA>  4.129551
26        <NA>  4.286294
27        <NA>  4.258179
28           Y  3.918690
29        <NA>  3.752468
30        <NA>  4.070887
31           Y  4.058540
32        <NA>  4.308313
33        <NA>  4.254858
34        <NA>  3.650220
35        <NA>  3.952759
36           Y  4.138911
37        <NA>  3.841253
38        <NA>  4.053278
39        <NA>  3.952263
40        <NA>  3.831942
41        <NA>  3.938618
42        <NA>  3.959048
43        <NA>  3.930177
44        <NA>  3.782907
45        <NA>  3.867163
46        <NA>  4.136693
47        <NA>  3.823215
48        <NA>  4.041099
49        <NA>  3.880098
50        <NA>  3.669454
51           Y  4.029818
52        <NA>  3.764323
53        <NA>  4.072553
54        <NA>  4.139595
55        <NA>  4.036563
56           Y  4.086740
57        <NA>  4.089358
58        <NA>  4.032313
59        <NA>  3.990262
60        <NA>  4.345966
61           Y  4.139036
62        <NA>  3.912880
63        <NA>  4.041223
64        <NA>  3.865279
65        <NA>  3.975427
66        <NA>  3.968791
67        <NA>  4.105275
68        <NA>  3.805072
69        <NA>  4.089852
70           Y  3.813947
71        <NA>  4.154348
72        <NA>  3.992997
73           Y  4.058475
74        <NA>  4.093328
75           Y  4.019205
76           Y  4.152910
77           Y  4.175624
78           Y  4.311996
79        <NA>  4.059912
80           Y  3.966214
81        <NA>  3.940277
82        <NA>  4.379176
83        <NA>  4.144147
84           Y  3.960076
85           Y  3.924314
86           Y  4.076846
87        <NA>  4.171999
88        <NA>  4.215812
89           Y  3.821951
90        <NA>  4.119819
91           Y  3.824250
92           Y  4.078508
93           Y  3.935337
94           Y  3.830297
95        <NA>  3.969295
96           Y  3.931549
97           Y  3.526715
98        <NA>  4.003508
99        <NA>  3.880599
100          Y  3.920924
101       <NA>  4.022677
102          Y  4.101613
103       <NA>  3.960888
104       <NA>  4.064412
105       <NA>  3.834159
106          Y  4.219119
107       <NA>  3.967249
108       <NA>  3.944910
109          Y  3.930995
110       <NA>  3.989005
111       <NA>  4.351425
112       <NA>  4.066976
113       <NA>  3.927269
114          Y  3.938739
115          Y  4.271694
116       <NA>  3.972534
117       <NA>  4.150076
118       <NA>  3.928875
119          Y  4.058829
120       <NA>  3.916164
121          Y  3.824933
122       <NA>  3.748669
123       <NA>  3.897229
124       <NA>  4.140008
125       <NA>  4.280480
126          Y  4.158507
127          Y  4.384277
128          Y  4.396899
129       <NA>  4.218378
130          Y  4.368716
131       <NA>  4.227045
132          Y  4.086559
133          Y  3.664094
134       <NA>  4.192900
135          Y  4.195769
136       <NA>  4.010934
137          Y  4.116780
138          Y  3.923291
139       <NA>  3.886227
140       <NA>  4.033894
141          Y  4.033137
142       <NA>  3.789624
143       <NA>  4.233798
144       <NA>  3.850588
145          Y  4.236659
146       <NA>  4.563047
147          Y  3.900137
148       <NA>  4.264080
149       <NA>  4.048912
150       <NA>  4.099126
151       <NA>  3.999838
152       <NA>  3.903822
153       <NA>  3.535851
154       <NA>  4.295383
155          Y  4.061669
156       <NA>        NA
157       <NA>        NA
158       <NA>        NA
159       <NA>        NA
160       <NA>        NA
161       <NA>        NA
162       <NA>        NA
163       <NA>        NA
164       <NA>        NA
165       <NA>        NA
166       <NA>        NA
167       <NA>        NA
168       <NA>        NA
169       <NA>        NA
170       <NA>        NA
171       <NA>        NA
172       <NA>        NA
173       <NA>        NA
174       <NA>        NA
175       <NA>        NA
176       <NA>        NA
177       <NA>        NA
178       <NA>        NA
179       <NA>        NA
180       <NA>        NA
181       <NA>        NA
182       <NA>        NA
183       <NA>        NA
184       <NA>        NA
185       <NA>        NA
186       <NA>        NA
187       <NA>        NA
188       <NA>        NA
189       <NA>        NA
190       <NA>        NA
191       <NA>        NA
192       <NA>        NA
193       <NA>        NA
194       <NA>        NA
195       <NA>        NA
196       <NA>        NA
197       <NA>        NA
198       <NA>        NA
199       <NA>        NA
200       <NA>        NA
201       <NA>        NA
202       <NA>        NA
203       <NA>        NA
204       <NA>        NA
205       <NA>        NA
206       <NA>        NA
207       <NA>        NA
208       <NA>        NA
209       <NA>        NA
210       <NA>        NA
211       <NA>        NA
212       <NA>        NA
213       <NA>        NA
214       <NA>        NA
215       <NA>        NA
216       <NA>  3.991518
217       <NA>  3.871037
218       <NA>  4.111620
219       <NA>  4.295654
220       <NA>  3.844464
221       <NA>  3.709591
222       <NA>  4.116721
223       <NA>  3.870917
224       <NA>  4.081825
225       <NA>  4.067721
226       <NA>  4.472735
227       <NA>  4.335865
228       <NA>  4.199849
229       <NA>  4.456463
230       <NA>  4.057840
231       <NA>  3.863950
232       <NA>  3.964976
233          N  3.903375
234       <NA>  4.155100
235       <NA>  4.149848
236       <NA>  3.972280
237       <NA>  4.434937
238       <NA>  4.012542
239       <NA>  4.556279
240       <NA>  4.545502
241       <NA>  4.189870
242       <NA>  4.507618
243       <NA>  4.176324
244       <NA>  4.111282
245       <NA>  4.428505
246       <NA>  4.097193
247       <NA>  3.950230
248       <NA>  4.317527
249       <NA>  4.182111
250       <NA>  4.396654
251       <NA>  3.839319
252       <NA>  4.310906
253       <NA>  4.078813
254       <NA>  4.380617
255       <NA>  3.856398
256       <NA>  4.590169
257       <NA>  4.179395
258       <NA>  4.035345
259       <NA>  3.930816
260       <NA>  4.212682
261       <NA>  4.155242
262       <NA>  3.993742
263       <NA>  4.288010
264       <NA>  4.059431
265       <NA>  4.338470
266       <NA>  3.956925
267       <NA>  4.328111
268       <NA>  4.349685
269       <NA>  4.333354
270       <NA>  3.865040
271       <NA>  4.130006
272       <NA>  4.145059
273       <NA>  4.237650
274          N  4.095722
275       <NA>  3.976744
276       <NA>  4.511975
277       <NA>  4.499879
278       <NA>  4.129845
279       <NA>  4.089548
280       <NA>  4.606000
281       <NA>  4.031684
282       <NA>  3.883730
283       <NA>  3.822597
284       <NA>  3.720674
285       <NA>  4.173343
286       <NA>  4.261594
287       <NA>  4.138022
288       <NA>  4.218968
289       <NA>  3.991661
290       <NA>  3.945646
291       <NA>  4.404076
292       <NA>  4.180484
293       <NA>  4.231135
294       <NA>  4.052863
295       <NA>  3.941177
296       <NA>  3.840258
297       <NA>  3.954041
298       <NA>  4.122124
299       <NA>  3.855313
300       <NA>  3.977862
301       <NA>  4.003031
302       <NA>  4.121179
303       <NA>  4.325789
304       <NA>  4.006242
305       <NA>  4.161272
306       <NA>  4.079270
307       <NA>  4.306855
308       <NA>  4.326332
309       <NA>  4.260949
310       <NA>  4.248080
311       <NA>  4.012514
312       <NA>  4.347608
313       <NA>  4.017278
314       <NA>  4.182197
315       <NA>  4.251114
316       <NA>  3.858479
317       <NA>  4.440506
318       <NA>  4.028412
319       <NA>  3.972269
320       <NA>  4.434554
321       <NA>  3.952204
322       <NA>  4.103837
323       <NA>  4.110850
324       <NA>  4.117688
325       <NA>  4.013520
326       <NA>  3.699630
327       <NA>  3.882692
328       <NA>  4.107256
329       <NA>  4.082395
330       <NA>  3.982675
331       <NA>  4.117749
332       <NA>  4.006744
333       <NA>  3.735111
334       <NA>  4.244293
335       <NA>  3.967445
336       <NA>  4.283030
337       <NA>  4.334178
338       <NA>  3.981520
339       <NA>  4.376559
340       <NA>  4.292818
341       <NA>  4.175456
342          2  4.266109
343       <NA>        NA
344       <NA>        NA
345       <NA>        NA
346       <NA>        NA
347       <NA>        NA
348       <NA>        NA
349       <NA>        NA
350       <NA>        NA
351       <NA>        NA
352       <NA>        NA
353       <NA>        NA
354       <NA>        NA
355       <NA>        NA
356       <NA>        NA
357       <NA>        NA
358       <NA>        NA
359       <NA>        NA
360       <NA>        NA
361       <NA>        NA
362       <NA>        NA
363       <NA>        NA
364       <NA>        NA
365       <NA>        NA
366       <NA>        NA
367       <NA>        NA
368       <NA>        NA
369       <NA>        NA
370       <NA>        NA
371       <NA>        NA
372       <NA>        NA
373       <NA>        NA
374       <NA>        NA
375       <NA>        NA
376       <NA>        NA
377       <NA>        NA
378       <NA>        NA
379       <NA>        NA
380       <NA>        NA
381       <NA>        NA
382       <NA>  4.814833
383       <NA>  4.474155
384       <NA>  4.515089
385       <NA>  4.583145
386       <NA>  4.578887
387       <NA>  4.640122
388       <NA>  4.672834
389       <NA>  4.589637
390       <NA>  4.451440
391       <NA>  4.366315
392       <NA>  4.476976
393       <NA>  4.280075
394       <NA>  4.493898
395       <NA>  4.631644
396       <NA>  4.575230
397       <NA>  4.327960
398       <NA>  4.504486
399       <NA>  4.678829
400       <NA>  4.373619
401       <NA>  4.566397
402       <NA>  4.390206
403       <NA>  4.119671
404       <NA>  4.531339
405       <NA>  4.182023
406       <NA>  4.634349
407       <NA>  4.328031
408       <NA>  4.224312
409       <NA>  4.521034
410       <NA>  4.566100
411       <NA>  4.243038
412       <NA>  4.552906
413       <NA>  4.325147
414       <NA>  4.321530
415       <NA>  4.264327
416       <NA>  4.626002
417       <NA>  4.225174
418       <NA>  4.219683
419       <NA>  4.583588
420       <NA>  4.517128
421       <NA>  4.357075
422       <NA>  4.425354
423       <NA>  4.323822
424       <NA>  4.500915
425       <NA>  4.695080
426       <NA>  4.572502
427       <NA>  4.835170
428       <NA>  4.541492
429       <NA>  4.424732
430       <NA>  4.158129
431       <NA>  4.795992
432       <NA>  4.390049
433       <NA>  4.410094
434       <NA>  4.309072
435       <NA>  4.620464
436       <NA>  4.367026
437       <NA>  4.443360
438       <NA>  4.613040
439       <NA>  4.441056
440       <NA>  4.745139
441       <NA>  4.245611
442       <NA>  4.292655
443       <NA>  4.710099
444       <NA>  4.421594
445       <NA>  4.067040
446       <NA>  4.853737
447       <NA>  4.253185
448       <NA>  4.448891
449       <NA>  4.212911
450       <NA>  4.494361
451       <NA>  4.402478
452       <NA>  4.199300
453       <NA>  4.789857
454       <NA>  4.211538
455       <NA>  4.508822
456       <NA>  4.605823
457       <NA>  4.692052
458       <NA>  4.537840
459       <NA>  4.360868
460       <NA>  4.396841
461       <NA>  4.464069
462       <NA>  4.765345
463       <NA>  4.731887
464       <NA>  4.317893
465       <NA>  4.620833
466       <NA>  4.591956
467       <NA>  4.995672
468       <NA>  5.003814
469       <NA>  4.823220
470       <NA>  4.605890
471       <NA>  4.578153
472       <NA>  4.885513
473       <NA>  4.906561
474       <NA>  4.914807
475       <NA>  4.734950
476       <NA>        NA
477       <NA>        NA
478       <NA>        NA
479       <NA>        NA
480       <NA>        NA
481       <NA>        NA
482       <NA>        NA
483       <NA>        NA
484       <NA>        NA
485       <NA>        NA
486       <NA>        NA
487       <NA>        NA
488       <NA>        NA
489       <NA>        NA
490       <NA>        NA
491       <NA>        NA
492       <NA>        NA
493       <NA>        NA
494       <NA>        NA
495       <NA>        NA
496       <NA>        NA
497       <NA>        NA
498       <NA>        NA
499       <NA>        NA
500       <NA>        NA
501       <NA>        NA
502       <NA>        NA
503       <NA>        NA
504       <NA>        NA
505       <NA>        NA
506       <NA>        NA
507       <NA>        NA
508       <NA>        NA
509       <NA>        NA
510       <NA>        NA
511       <NA>        NA
512       <NA>        NA
513       <NA>        NA
514       <NA>        NA
515       <NA>  5.700533
516       <NA>  5.600874
517       <NA>  5.640983
518       <NA>  5.349059
519       <NA>  5.349432
520       <NA>  5.157335
521       <NA>  5.480632
522       <NA>  5.568130
523       <NA>  5.598492
524       <NA>  5.632051
525       <NA>  5.352134
526       <NA>  5.663278
527       <NA>  5.451533
528       <NA>  5.157956
529       <NA>  5.273777
530       <NA>  5.222730
531       <NA>  5.189449
532       <NA>  5.171091
533       <NA>  5.194670
534       <NA>  5.434768
535       <NA>  5.232679
536       <NA>  5.125214
537       <NA>  5.079366
538       <NA>  5.041921
539       <NA>  5.082263
540       <NA>  5.415455
541       <NA>  5.157884
542       <NA>  5.252326
543       <NA>  5.432417
544       <NA>  5.239270
545       <NA>  5.275286
546       <NA>  5.155576
547       <NA>  5.405964
548       <NA>  5.534177
549       <NA>  5.259326
550       <NA>  5.073481
551       <NA>  5.418345
552       <NA>  5.471626
553       <NA>  5.195893
554       <NA>  5.538252
555       <NA>  5.459515
556       <NA>  5.468355
557       <NA>  5.380998
558       <NA>  5.065075
559       <NA>  5.462305
560       <NA>  5.402515
561       <NA>  5.146772
562       <NA>  5.040351
563       <NA>  5.407088
564       <NA>  5.233550
565       <NA>  5.747026
566       <NA>  5.233359
567       <NA>  5.648286
568       <NA>  5.464880
569       <NA>  5.096133
570       <NA>  5.277781
571       <NA>  5.539862
572       <NA>  5.159951
573       <NA>  5.195746
574       <NA>  5.237172
575       <NA>  5.180285
576       <NA>  5.180749
577       <NA>  5.764127
578       <NA>  5.154133
579       <NA>  5.483589
580       <NA>  5.401604
581       <NA>  5.488415
582       <NA>  5.173245
583       <NA>  5.129888
584       <NA>  6.215643
585       <NA>  6.134307
586       <NA>  6.044461
587       <NA>  6.212898
588       <NA>  6.525253
589       <NA>  6.469152
590       <NA>  6.327194
591       <NA>  6.291526
592       <NA>  6.230008
593       <NA>        NA
594       <NA>        NA
595       <NA>        NA
596       <NA>        NA
597       <NA>        NA
598       <NA>        NA
599       <NA>        NA
600       <NA>        NA
601       <NA>        NA
602       <NA>        NA
603       <NA>        NA
604       <NA>        NA
605       <NA>        NA
606       <NA>        NA
607       <NA>        NA
608       <NA>        NA
609       <NA>        NA
610       <NA>        NA
611          N  6.756172
612       <NA>  6.519212
613       <NA>  6.325626
614       <NA>  6.245649
615       <NA>  6.565046
616          Y  6.591842
617       <NA>  6.320392
618       <NA>  6.313265
619          N  6.520100
620          N  6.461723
621       <NA>  6.389255
622       <NA>  6.376899
623       <NA>  5.874616
624       <NA>  6.079242
625       <NA>  5.985972
626       <NA>  6.077646
627       <NA>  6.296338
628       <NA>  5.861135
629          Y  6.223682
630       <NA>  6.229736
631       <NA>  6.110690
632       <NA>  6.236400
633       <NA>  6.098210
634       <NA>  6.411226
635       <NA>  6.439952
636       <NA>  5.687091
637       <NA>  6.256476
638       <NA>  6.142738
639       <NA>  5.937312
640          Y  6.185021
641          Y  6.296402
642       <NA>  6.048440
643          Y  6.104892
644       <NA>  6.153833
645       <NA>  6.192625
646       <NA>  6.154033
647       <NA>  5.836078
648       <NA>  5.874363
649       <NA>  6.027413
650       <NA>  6.477395
651       <NA>  6.009503
652       <NA>  5.873449
653       <NA>  5.922556
654       <NA>  6.310396
655       <NA>  5.983001
656       <NA>  6.036799
657       <NA>  6.324921
658       <NA>  6.281023
659       <NA>  6.049633
660       <NA>  5.958667
661       <NA>  6.200817
662       <NA>  6.272608
663       <NA>  5.953289
664          Y  6.259054
665       <NA>  6.174896
666       <NA>  6.138049
667       <NA>  7.649965
668       <NA>  7.532024
669       <NA>  7.425346
670       <NA>  7.549838
671       <NA>  7.451369
672       <NA>  7.513671
673       <NA>  7.681176
674          N  7.600269
675       <NA>        NA
676       <NA>        NA
677       <NA>        NA
678       <NA>        NA
679       <NA>        NA
680       <NA>        NA
681       <NA>        NA
682       <NA>        NA
683       <NA>        NA
684       <NA>        NA
685       <NA>        NA
686       <NA>        NA
687       <NA>        NA
688       <NA>        NA
689       <NA>        NA
690       <NA>        NA
691       <NA>        NA
692       <NA>        NA
693       <NA>        NA
694       <NA>        NA
695       <NA>        NA
696       <NA>        NA
697       <NA>        NA
698       <NA>        NA
699       <NA>        NA
700       <NA>        NA
701       <NA>        NA
702       <NA>  8.538176
703       <NA>  8.600173
704       <NA>  8.439642
705       <NA>  8.813960
706       <NA>  8.538191
707       <NA>  8.168513
708       <NA>  7.932521
709       <NA>  8.109984
710       <NA>  7.858755
711       <NA>  8.113696
712       <NA>  7.872023
713       <NA>  7.829224
714       <NA>  7.736618
715       <NA>  8.441267
716       <NA>  8.034820
717       <NA>  7.590457
718       <NA>  8.122827
719       <NA>  8.110128
720       <NA>  7.968968
721       <NA>  7.801010
722       <NA>  7.849104
723       <NA>  7.442824
724       <NA>  8.166012
725       <NA>  7.864415
726       <NA>  8.177747
727       <NA>  8.008055
728       <NA>  7.902348
729       <NA>  8.057404
730       <NA>  7.889184
731       <NA>  7.901946
732       <NA>  8.213785
733       <NA>  7.859453
734       <NA>  8.293056
735       <NA>  7.903874
736       <NA>  7.757460
737       <NA>  7.977984
738       <NA>  7.873699
739       <NA>  8.099806
740       <NA>  7.959608
741       <NA>  8.005570
742       <NA>  7.991690
743       <NA>  7.629156
744       <NA>  7.930108
745       <NA>  8.154555
746       <NA>  7.893755
747       <NA>  8.045484
748       <NA>  7.856610
749       <NA>  7.683940
750       <NA> 11.360772
751       <NA> 11.100219
752       <NA> 10.862076
753          Y 11.016912
754       <NA> 11.056551
755       <NA>        NA
756       <NA>        NA
757       <NA>        NA
758       <NA>        NA
759       <NA>        NA
760       <NA>        NA
761       <NA>        NA
762       <NA>        NA
763       <NA>        NA
764       <NA>        NA
765       <NA>        NA
766       <NA>        NA
767       <NA>        NA
768       <NA>        NA
769       <NA>        NA
770       <NA>        NA
771       <NA>        NA
772       <NA>        NA
773       <NA> 10.910503
774       <NA> 10.952710
775       <NA> 11.010630
776       <NA> 10.832632
777       <NA> 10.804228
778       <NA>  9.749238
779       <NA>  9.942216
780       <NA> 10.127860
781       <NA>  9.798381
782          Y  9.789847
783       <NA>  9.614897
784       <NA>  9.728010
785       <NA>  9.900402
786       <NA>  9.895083
787       <NA>  9.927820
788       <NA>  9.916280
789       <NA>  9.847538
790       <NA>  9.873366
791       <NA> 10.145782
792       <NA>  9.927355
793       <NA>  9.905220
794       <NA>  9.953278
795          Y  9.861769
796       <NA>  9.838720
797       <NA>  9.950482
798       <NA>  9.882120
799       <NA>  9.916330
800       <NA> 10.076311
801       <NA> 10.059305
802       <NA>  9.965447
803       <NA>  9.804127
804       <NA>  9.751070
805       <NA>  9.792538
806       <NA>  9.871957
807       <NA>  9.994643
808       <NA> 10.070241
809       <NA>  9.851530
810       <NA>  9.654886
811       <NA> 10.066954
812       <NA>  9.837881
813       <NA>  9.863769
814       <NA>  9.979368
815       <NA> 10.213945
816       <NA> 10.115429
817       <NA> 14.416301
818       <NA> 14.197060
819       <NA> 14.669467
820       <NA> 14.664155
821       <NA>        NA
822       <NA>        NA
823       <NA>        NA
824       <NA>        NA
825       <NA>        NA
826       <NA>        NA
827       <NA>        NA
828       <NA>        NA
829       <NA>        NA
830       <NA>        NA
831       <NA>        NA
832       <NA>        NA
833       <NA>        NA
834       <NA>        NA
835       <NA>        NA
836       <NA>        NA
837       <NA>        NA
838       <NA>        NA
839       <NA> 14.157665
840       <NA> 14.026238
841       <NA> 13.798979
842       <NA> 13.980592
843       <NA> 12.528002
844       <NA> 12.507252
845       <NA> 12.583383
846       <NA> 12.626965
847       <NA> 12.540796
848       <NA> 12.313872
849       <NA> 12.955107
850       <NA> 12.503455
851       <NA> 12.425007
852       <NA> 12.718740
853       <NA> 12.709809
854       <NA> 12.554592
855       <NA> 12.435392
856       <NA> 12.496626
857       <NA> 12.320248
858       <NA> 12.393628
859       <NA> 12.679804
860       <NA> 12.175063
861       <NA> 12.563563
862       <NA> 12.680667
863       <NA> 12.779860
864       <NA> 12.373359
865       <NA> 12.706026
866       <NA> 12.805650
867       <NA> 12.625226
868       <NA> 12.725137
869       <NA> 12.619746
870       <NA> 12.485814
871       <NA> 12.438928
872       <NA> 12.407697
873       <NA> 12.444207
874       <NA> 12.747430
875       <NA> 12.553742
876       <NA> 12.439180
877       <NA> 12.716731
878       <NA> 12.468676
879       <NA> 12.648160
880       <NA> 19.546396
881       <NA> 19.095462
882       <NA> 19.124670
883       <NA> 19.454349
884       <NA>        NA
885       <NA>        NA
886       <NA>        NA
887       <NA>        NA
888       <NA>        NA
889       <NA>        NA
890       <NA>        NA
891       <NA>        NA
892       <NA>        NA
893       <NA>        NA
894       <NA>        NA
895       <NA>        NA
896       <NA>        NA
897       <NA>        NA
898       <NA>        NA
899       <NA>        NA
900       <NA>        NA
901       <NA>        NA
902       <NA> 18.244533
903       <NA> 17.707753
904       <NA> 17.872669
905       <NA> 15.818165
906       <NA> 15.745895
907          N 16.019507
908       <NA> 16.103556
909       <NA> 16.202645
910       <NA> 15.934590
911       <NA> 16.215263
912       <NA> 15.742315
913       <NA> 15.737214
914       <NA> 16.161898
915       <NA> 16.195697
916       <NA> 16.060271
917       <NA> 16.115513
918       <NA> 15.641730
919       <NA> 15.903854
920       <NA> 15.869049
921       <NA> 16.087205
922       <NA> 15.928673
923       <NA> 16.096859
924       <NA> 15.912562
925       <NA> 15.990284
926       <NA> 16.034631
927       <NA> 16.167671
928       <NA> 15.919300
929       <NA> 15.872934
930       <NA> 16.336027
931       <NA> 15.962674
932       <NA> 16.058532
933       <NA> 16.172012
934       <NA> 16.404351
935       <NA> 16.218213
936       <NA> 15.807350
937          N 15.939829
938       <NA> 16.026186
939       <NA> 15.905642
940       <NA> 15.756689
941       <NA> 25.447177
942       <NA> 25.349509
943       <NA> 25.306073
944       <NA> 25.376368
945       <NA>        NA
946       <NA>        NA
947       <NA>        NA
948       <NA>        NA
949       <NA>        NA
950       <NA>        NA
951       <NA>        NA
952       <NA>        NA
953       <NA>        NA
954       <NA>        NA
955       <NA>        NA
956       <NA>        NA
957       <NA>        NA
958       <NA>        NA
959       <NA>        NA
960       <NA>        NA
961       <NA>        NA
962       <NA>        NA
963       <NA> 22.989226
964       <NA> 23.010487
965       <NA> 22.974484
966       <NA> 20.256620
967       <NA> 20.269787
968       <NA> 20.629028
969       <NA> 20.489896
970       <NA> 20.374301
971          Y 20.159805
972       <NA> 20.461975
973       <NA> 20.326244
974       <NA> 20.221890
975       <NA> 20.234575
976       <NA> 20.355661
977       <NA> 20.172728
978       <NA> 20.458419
979          Y 20.360573
980       <NA> 20.317379
981       <NA> 20.159305
982       <NA> 20.087718
983       <NA> 20.110715
984       <NA> 19.998388
985       <NA> 20.108398
986       <NA> 20.213543
987       <NA> 20.122793
988       <NA> 20.669172
989       <NA> 20.331733
990       <NA> 20.161959
991       <NA> 20.248634
992       <NA> 20.162556
993       <NA> 20.519142
994       <NA> 20.259703
995       <NA> 19.970717
996       <NA> 20.317138
997       <NA> 20.385787
998       <NA> 20.245450
999       <NA> 20.204008
1000      <NA> 20.493939
1001      <NA> 20.411062
1002      <NA> 32.965480
1003      <NA> 32.721821
1004      <NA> 33.240541
1005      <NA>        NA
1006      <NA>        NA
1007      <NA>        NA
1008      <NA>        NA
1009      <NA>        NA
1010      <NA>        NA
1011      <NA>        NA
1012      <NA>        NA
1013      <NA>        NA
1014      <NA>        NA
1015      <NA>        NA
1016      <NA>        NA
1017      <NA>        NA
1018      <NA>        NA
1019      <NA>        NA
1020      <NA>        NA
1021      <NA>        NA
1022      <NA>        NA
1023      <NA> 29.651782
1024      <NA> 25.320905
1025      <NA> 25.277072
1026      <NA> 25.426995
1027      <NA> 25.649796
1028      <NA> 25.443149
1029      <NA> 25.723984
1030      <NA> 25.390111
1031      <NA> 25.413160
1032      <NA> 25.251020
1033      <NA> 25.243851
1034      <NA> 25.193123
1035      <NA> 25.571238
1036      <NA> 25.424359
1037      <NA> 25.583007
1038      <NA> 25.401022
1039      <NA> 25.509987
1040      <NA> 25.324085
1041      <NA> 25.248657
1042      <NA> 25.512930
1043      <NA> 25.270236
1044      <NA> 25.347544
1045      <NA> 25.185450
1046      <NA> 25.198018
1047      <NA> 25.403791
1048      <NA> 25.357156
1049      <NA> 25.296619
1050      <NA> 25.532701
1051      <NA> 25.242244
1052      <NA> 25.104720
1053      <NA> 25.546274
1054      <NA> 25.184024
1055      <NA> 25.548208
1056      <NA> 25.535378
1057      <NA> 25.209085
1058      <NA> 25.684138
1059      <NA> 25.328171
1060      <NA> 42.288193
1061      <NA>        NA
1062      <NA>        NA
1063      <NA>        NA
1064      <NA>        NA
1065      <NA>        NA
1066      <NA>        NA
1067      <NA>        NA
1068      <NA>        NA
1069      <NA>        NA
1070      <NA>        NA
1071      <NA>        NA
1072      <NA>        NA
1073      <NA>        NA
1074      <NA>        NA
1075      <NA>        NA
1076      <NA>        NA
1077      <NA>        NA
1078      <NA>        NA
1079      <NA>        NA
1080      <NA>        NA
1081      <NA>        NA
1082      <NA>        NA
1083      <NA>        NA
1084      <NA>        NA
1085      <NA>        NA
1086      <NA>        NA
1087      <NA>        NA
1088      <NA>        NA
1089      <NA>        NA
1090      <NA>        NA
1091      <NA>        NA
1092      <NA>        NA
1093      <NA>        NA
1094      <NA>        NA
1095      <NA>        NA
1096      <NA>        NA
1097      <NA>        NA
1098      <NA>        NA
1099      <NA>        NA
1100      <NA>        NA
1101      <NA>        NA
1102      <NA>        NA
1103      <NA>        NA
1104      <NA>        NA
1105      <NA>        NA
1106      <NA>        NA
1107      <NA>        NA
1108      <NA>        NA
1109      <NA>        NA
1110      <NA>        NA
1111      <NA>        NA
1112      <NA>        NA
1113      <NA>        NA
1114      <NA>        NA
1115      <NA>        NA
1116      <NA>        NA
1117      <NA>        NA
1118      <NA>        NA
1119      <NA>        NA
1120      <NA>        NA
1121      <NA>        NA
1122      <NA>        NA
1123      <NA>        NA
1124      <NA>        NA
1125      <NA>        NA
1126      <NA>        NA
1127      <NA>        NA
1128      <NA>        NA
1129      <NA>        NA
1130      <NA>        NA
1131      <NA>        NA
1132      <NA>        NA
1133      <NA>        NA
1134      <NA>        NA
1135      <NA>        NA
1136      <NA>        NA
1137      <NA>        NA
1138      <NA>        NA
1139      <NA>        NA
1140      <NA>        NA
1141      <NA>        NA
1142      <NA>        NA
1143      <NA> 50.001149
1144      <NA> 49.831917
1145      <NA> 49.812766
1146      <NA> 50.092470
1147      <NA> 49.858388
1148      <NA> 49.968799
1149      <NA> 49.891059
1150      <NA> 50.176740
1151      <NA> 49.907082
1152      <NA> 49.706507
1153      <NA> 50.184944
1154      <NA> 50.126574
1155      <NA> 49.462248
1156      <NA> 49.920601
1157      <NA> 49.899879
1158      <NA> 49.718128
1159      <NA> 50.066502
1160      <NA> 49.848419
1161      <NA> 50.009365
1162      <NA> 49.919848
1163      <NA> 49.694198
1164      <NA> 49.532585
1165      <NA> 50.106733
1166      <NA> 49.896587
1167      <NA> 49.760386
1168      <NA> 49.920258
1169      <NA> 50.018159
1170      <NA> 49.770226
1171      <NA> 50.001088
1172      <NA> 49.687627
1173      <NA>        NA
1174      <NA>        NA
1175      <NA>        NA
1176      <NA>        NA
1177      <NA>        NA
1178      <NA>        NA
1179      <NA>        NA
1180      <NA>        NA
1181      <NA>        NA
1182      <NA>        NA
1183      <NA>        NA
1184      <NA>        NA
1185      <NA>        NA
1186      <NA>        NA
1187      <NA>        NA
1188      <NA>        NA
1189      <NA>        NA
1190      <NA>        NA
1191      <NA>        NA
1192      <NA>        NA
1193      <NA>        NA
1194      <NA>        NA
1195      <NA>        NA
1196      <NA>        NA
1197      <NA>        NA
1198      <NA>        NA
1199      <NA>        NA
1200      <NA>        NA
1201      <NA>        NA
1202      <NA>        NA
1203      <NA>        NA
1204      <NA>        NA
1205      <NA>        NA
1206      <NA>        NA
1207      <NA>        NA
1208      <NA>        NA
1209      <NA>        NA
1210      <NA>        NA
1211      <NA>        NA
1212      <NA>        NA
1213      <NA>        NA
1214      <NA>        NA
1215      <NA>        NA
1216      <NA>        NA
1217      <NA>        NA
1218      <NA>        NA
1219      <NA>        NA
1220      <NA>        NA
1221      <NA>        NA
1222      <NA>        NA
1223      <NA>        NA
1224      <NA>        NA
1225      <NA>        NA
1226      <NA>        NA
1227      <NA>        NA
1228      <NA>        NA
1229      <NA>        NA
1230      <NA>        NA
1231      <NA>        NA
1232      <NA>        NA
1233      <NA>        NA
1234      <NA>        NA
1235      <NA>        NA
1236      <NA>        NA
1237      <NA>        NA
1238      <NA>        NA
1239      <NA>        NA
1240      <NA>        NA
1241      <NA>        NA
1242      <NA>        NA
1243      <NA>        NA
1244      <NA>        NA
1245      <NA>        NA
1246      <NA>        NA
1247      <NA>        NA
1248      <NA>        NA
1249      <NA>        NA
1250      <NA>        NA
1251      <NA>        NA
1252      <NA>        NA
1253      <NA>        NA
1254      <NA>        NA
1255      <NA>        NA
1256      <NA> 63.528494
1257      <NA> 63.760973
1258      <NA> 63.696884
1259      <NA> 63.632992
1260      <NA> 63.376237
1261      <NA> 64.065534
1262      <NA> 63.857197
1263      <NA> 63.562629
1264      <NA> 63.495964
1265      <NA> 63.692196
1266      <NA> 63.590003
1267      <NA> 63.350309
1268      <NA> 63.690672
1269      <NA> 63.520392
1270      <NA> 63.692243
1271      <NA> 63.502841
1272      <NA> 63.609270
1273      <NA> 63.752470
1274      <NA> 63.492699
1275      <NA> 63.745796
1276      <NA> 63.543757
1277      <NA> 63.615452
1278      <NA> 63.753489
1279      <NA> 63.884817
1280      <NA> 63.926851
1281      <NA> 63.521275
1282      <NA> 63.418812
1283      <NA> 63.764012
1284      <NA>        NA
1285      <NA>        NA
1286      <NA>        NA
1287      <NA>        NA
1288      <NA>        NA
1289      <NA>        NA
1290      <NA>        NA
1291      <NA>        NA
1292      <NA>        NA
1293      <NA>        NA
1294      <NA>        NA
1295      <NA>        NA
1296      <NA>        NA
1297      <NA>        NA
1298      <NA>        NA
1299      <NA>        NA
1300      <NA>        NA
1301      <NA>        NA
1302      <NA>        NA
1303      <NA>        NA
1304      <NA> 64.434747
1305      <NA> 64.975639
1306      <NA> 64.812816
1307      <NA> 64.822377
1308      <NA> 64.644733
1309      <NA> 64.340282
1310      <NA> 64.473237
1311      <NA> 64.699403
1312      <NA> 65.004615
1313      <NA> 64.722035
1314      <NA> 64.386538
1315      <NA> 64.667017
1316      <NA> 64.907174
1317      <NA> 64.520354
1318      <NA> 64.686133
1319      <NA> 64.894768
1320      <NA> 64.657335
1321      <NA> 64.773103
1322      <NA> 64.524492
1323      <NA> 64.331944
1324      <NA> 64.803224
1325      <NA> 64.358045
1326      <NA>        NA
1327      <NA>        NA
1328      <NA>        NA
1329      <NA>        NA
1330      <NA>        NA
1331      <NA>        NA
1332      <NA>        NA
1333      <NA>        NA
1334      <NA>        NA
1335      <NA>        NA
1336      <NA>        NA
1337      <NA>        NA
1338      <NA>        NA
1339      <NA>        NA
1340      <NA>        NA
1341      <NA>        NA
1342      <NA>        NA
1343      <NA>        NA
1344      <NA>        NA
1345      <NA>        NA
1346      <NA> 64.985001
1347      <NA> 65.321814
1348      <NA> 65.281855
1349      <NA> 65.475168
1350      <NA> 65.173232
1351      <NA> 65.144291
1352      <NA> 65.180670
1353      <NA> 65.131083
1354      <NA> 65.396993
1355      <NA> 65.065100
1356      <NA> 65.014401
1357      <NA> 65.062862
1358      <NA> 65.126453
1359      <NA> 65.005888
1360      <NA> 65.165616
1361      <NA> 64.964095
1362      <NA> 65.186881
1363      <NA> 65.418700
1364      <NA> 64.747750
1365      <NA> 65.239854
1366      <NA>        NA
1367      <NA>        NA
1368      <NA>        NA
1369      <NA>        NA
1370      <NA>        NA
1371      <NA>        NA
1372      <NA>        NA
1373      <NA>        NA
1374      <NA>        NA
1375      <NA>        NA
1376      <NA>        NA
1377      <NA>        NA
1378      <NA>        NA
1379      <NA>        NA
1380      <NA>        NA
1381      <NA>        NA
1382      <NA>        NA
1383      <NA>        NA
1384      <NA>        NA
1385      <NA>        NA
1386      <NA> 65.602948
1387      <NA> 65.472074
1388      <NA> 65.377037
1389      <NA> 65.581256
1390      <NA> 65.554622
1391      <NA> 65.655166
1392      <NA> 65.305949
1393      <NA> 65.216083
1394      <NA> 65.271489
1395      <NA> 65.430919
1396      <NA> 65.447314
1397      <NA> 65.731145
1398      <NA> 65.155773
1399      <NA> 65.150051
1400      <NA> 65.453374
1401      <NA> 65.715952
1402      <NA> 65.312677
1403      <NA> 65.198242
1404      <NA> 64.947175
1405      <NA> 65.324938
1406      <NA>        NA

plot full predictions

  • What do you think if the predictions? Any better than model 1?
  • Do you agree with the Hmax values generated by the model?
  • If not, why do you think the model is generating those values? (also see next slide)
fit2.predictions %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop), alpha=.2) +
  geom_line(aes(y=prediction, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2", guide="none") +
  facet_wrap(~pop)

DPR data + all predictions

  • This may be helpful in answering the last question
fit2.predictions %>% ggplot(aes(x=day)) +
  geom_line(aes(group=plantID, y=height_cm, color=pop),
            alpha=.2,
            data = {fit2.predictions %>% filter(pop=="DPR")}) +
  geom_line(aes(y=prediction, color=pop)) +
  scale_color_brewer(type="qual", palette = "Set2") 

Let’s try a new model

  • Having Hmax alone vary by population did not produce a very good fit.
  • Exercise: Create and fit a model that allows Hmax and k parameters to vary by population
  • Plot your results: Any better?
  • Perform a model comparison among all 3 models using the loo technique shown above: which fits best?

My formula

f3 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + delta ~ 1,
                  Hmax + k ~ pop,
                  nl=TRUE)

get the priors

Please review your notes and figure out how to get the possible priors for the formula on the previous slide

                prior class      coef group resp dpar nlpar lb ub       source
 student_t(3, 0, 4.2) sigma                                  0         default
               (flat)     b                           delta            default
               (flat)     b Intercept                 delta       (vectorized)
               (flat)     b                            Hmax            default
               (flat)     b Intercept                  Hmax       (vectorized)
               (flat)     b     popCC                  Hmax       (vectorized)
               (flat)     b    popDPR                  Hmax       (vectorized)
               (flat)     b                            Hmin            default
               (flat)     b Intercept                  Hmin       (vectorized)
               (flat)     b                               k            default
               (flat)     b Intercept                     k       (vectorized)
               (flat)     b     popCC                     k       (vectorized)
               (flat)     b    popDPR                     k       (vectorized)

Set the priors

What priors should we use?

prior3 <- c(set_prior("normal(5,3)", lb=0, nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("normal(0, 1)", nlpar = "k", coef="popCC"),
            set_prior("normal(0, 1)", nlpar = "k", coef= "popDPR"),
            set_prior("gamma(20, 3)", nlpar="delta")

)

I use a gamma on k to keep the Intercept k positive.

I use a zero-centered normal for the popCC and popDPR k priors so that these can vary above or below the intercept value of k

model summary

Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         delta ~ 1
         Hmax ~ pop
         k ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      4.02      0.53     2.96     5.03 1.00     3116     2007
delta_Intercept     6.69      0.51     5.71     7.73 1.00     2418     2610
Hmax_Intercept     64.88      1.33    62.35    67.61 1.00     2464     2303
Hmax_popCC        -15.19     20.51   -44.82    27.93 1.00     2039     2773
Hmax_popDPR       -29.09      6.70   -38.48   -12.02 1.00     2102     1281
k_Intercept         0.67      0.01     0.65     0.69 1.00     2273     2254
k_popCC             0.11      0.10    -0.03     0.34 1.00     1918     2441
k_popDPR            0.32      0.07     0.18     0.44 1.00     2206     1281

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.17      0.26    10.68    11.70 1.00     4326     2927

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

model plot

pairs

You can use a pairs plot to see how to parameters co-vary during the fit. Generally we don’t want to see a strong correlation among parameters. Here we see that Hmax and K for popCC have an inverse correlation. This isn’t ideal but I think it is because we have so little CC data. We could impose narrower priors on these parameters but given the lack of data I think I will just leave it.

Predictions

Model comparison

fit1 <- add_criterion(fit1, "loo")
fit2 <- add_criterion(fit2, "loo")
fit3 <- add_criterion(fit3, "loo")
loo_compare(fit1, fit2, fit3)
     elpd_diff se_diff
fit3   0.0       0.0  
fit2  -8.0       3.9  
fit1 -18.5       6.9  

Next

For Monday March 4th

  • Fit a model where delta and Hmax can vary by pop.
  • Fit a model where delta, k, and Hmax can vary by pop. Don’t worry if this doesn’t work well (but if it doesn’t, think about why it doesn’t)
  • Do model diagnostics, plots, and predictive plots on the above.
  • Do a model comparison between all 5 models

Other stuff that Julin is thinking about for future meetings

  • Uncertainty on the predictions
  • On to random effects and all pops
  • CENSORED DATA?

Delta and Hmax by pop:

f4 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin + k ~ 1,
                  Hmax + delta ~ pop,
                  nl=TRUE)

Set the priors

What priors should we use?

prior4 <- c(set_prior("normal(5,3)", lb=0, nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("gamma(20, 3)", nlpar="delta"),
            set_prior("normal(0, 5)", nlpar = "delta", coef="popCC"),
            set_prior("normal(0, 5)", nlpar = "delta", coef= "popDPR")

)

model summary

Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         k ~ 1
         Hmax ~ pop
         delta ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      3.87      0.55     2.76     4.93 1.00     4015     2324
k_Intercept         0.67      0.01     0.65     0.69 1.00     2863     2707
Hmax_Intercept     64.89      1.31    62.36    67.53 1.00     3346     2812
Hmax_popCC         -2.41     13.93   -27.61    26.85 1.00     3004     2632
Hmax_popDPR        13.92     11.21    -6.90    37.17 1.00     3676     2696
delta_Intercept     6.76      0.50     5.84     7.82 1.00     2351     2695
delta_popCC        -1.87      1.13    -3.87     0.58 1.00     2642     2781
delta_popDPR       -3.21      0.67    -4.56    -1.92 1.00     2734     2594

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.17      0.27    10.66    11.71 1.00     5526     3014

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

model plot

Predictions

Model comparison

fit4 <- add_criterion(fit4, "loo")
loo_compare(fit1, fit2, fit3, fit4)
     elpd_diff se_diff
fit3   0.0       0.0  
fit4  -0.1       1.3  
fit2  -8.0       3.9  
fit1 -18.5       6.9  

Hmax, delta, and k

f5 <- brmsformula(height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta),
                  Hmin ~ 1,
                  Hmax + k + delta ~ pop,
                  nl=TRUE)

Set the priors

What priors should we use?

prior5 <- c(set_prior("normal(5,3)", lb=0, nlpar="Hmin"),
            set_prior("normal(60,15)", nlpar="Hmax"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popCC"),
            set_prior("normal(0,20)", nlpar="Hmax", coef="popDPR"),
            set_prior("gamma(2,2)", nlpar = "k"), # gamma constrains to be > 0
            set_prior("normal(0, 1)", nlpar = "k", coef="popCC"),
            set_prior("normal(0, 1)", nlpar = "k", coef= "popDPR"),
            set_prior("gamma(20, 3)", nlpar="delta"),
            set_prior("normal(0, 5)", nlpar = "delta", coef="popCC"),
            set_prior("normal(0, 5)", nlpar = "delta", coef= "popDPR")

)

model summary

Running /Library/Frameworks/R.framework/Resources/bin/R CMD SHLIB foo.c
using C compiler: ‘Apple clang version 13.1.6 (clang-1316.0.21.2.5)’
using SDK: ‘MacOSX12.3.sdk’
clang -arch arm64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/Rcpp/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/unsupported"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/BH/include" -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/src/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppParallel/include/"  -I"/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/rstan/include" -DEIGEN_NO_DEBUG  -DBOOST_DISABLE_ASSERTS  -DBOOST_PENDING_INTEGER_LOG2_HPP  -DSTAN_THREADS  -DUSE_STANC3 -DSTRICT_R_HEADERS  -DBOOST_PHOENIX_NO_VARIADIC_EXPRESSION  -D_HAS_AUTO_PTR_ETC=0  -include '/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp'  -D_REENTRANT -DRCPP_PARALLEL_USE_TBB=1   "-I/opt/R/arm64/include -I/opt/homebrew/include"    -fPIC  -falign-functions=64 -Wall -g -O2  -c foo.c -o foo.o
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:88:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:1: error: unknown type name 'namespace'
namespace Eigen {
^
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/src/Core/util/Macros.h:628:16: error: expected ';' after top level declarator
namespace Eigen {
               ^
               ;
In file included from <built-in>:1:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/StanHeaders/include/stan/math/prim/fun/Eigen.hpp:22:
In file included from /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Dense:1:
/Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library/RcppEigen/include/Eigen/Core:96:10: fatal error: 'complex' file not found
#include <complex>
         ^~~~~~~~~
3 errors generated.
make: *** [foo.o] Error 1
 Family: gaussian 
  Links: mu = identity; sigma = identity 
Formula: height_cm ~ Hmax - (Hmax - Hmin) * exp(-(k/100 * day)^delta) 
         Hmin ~ 1
         Hmax ~ pop
         k ~ pop
         delta ~ pop
   Data: datbcd (Number of observations: 907) 
  Draws: 4 chains, each with iter = 2000; warmup = 1000; thin = 1;
         total post-warmup draws = 4000

Population-Level Effects: 
                Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
Hmin_Intercept      3.93      0.55     2.81     4.99 1.00     3053     2664
Hmax_Intercept     64.74      1.30    62.25    67.30 1.00     3268     2751
Hmax_popCC         -7.90     21.50   -44.40    33.56 1.00     1161     2015
Hmax_popDPR       -17.32     18.16   -37.95    26.67 1.01      699     1101
k_Intercept         0.67      0.01     0.65     0.69 1.00     2759     2670
k_popCC             0.02      0.14    -0.25     0.32 1.00     1334     2132
k_popDPR            0.23      0.15    -0.06     0.45 1.01      724     1290
delta_Intercept     6.82      0.51     5.88     7.86 1.00     2605     2515
delta_popCC        -1.72      2.17    -4.61     4.25 1.00     2168     1430
delta_popDPR       -0.56      3.07    -4.19     7.27 1.00      860     1391

Family Specific Parameters: 
      Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sigma    11.17      0.26    10.68    11.70 1.00     3617     2462

Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
and Tail_ESS are effective sample size measures, and Rhat is the potential
scale reduction factor on split chains (at convergence, Rhat = 1).

model plot

DPR Pairs

Intercept Pairs

Predictions

Model comparison

fit5 <- add_criterion(fit5, "loo", moment_match=TRUE)
loo_compare(fit1, fit2, fit3, fit4, fit5)
     elpd_diff se_diff
fit3   0.0       0.0  
fit5  -0.1       0.9  
fit4  -0.1       1.3  
fit2  -8.0       3.9  
fit1 -18.5       6.9